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Abstract 

I examine the solution of the BFKL equation with NLO corrections relevant for deep inelastic 

scattering. Particular emphasis is placed on the part played by the running of the coupling. It is 

shown that the solution factorizes into a part describing the evolution in Q 2 , and a constant part 

describing the input distribution. The latter is infrared dominated, being described by a coupling 

which grows as x decreases, and thus being contaminated by infrared renormalons. Hence, for this 

part we agree with previous assertions that predictive power breaks down for small enough x at any 

Q 2 . However, the former is ultraviolet dominated, being described by a coupling which falls like 

1 

l/(ln(Q 2 /A 2 ) + A(a s {Q 2 ) ln(l/x)) 2 ) with decreasing x, and thus is perturbatively calculable at all 
x. Therefore, although the BFKL equation is unable to predict the input for a structure function for 
small x, it is able to predict its evolution in Q 2 , as we would expect from the factorization theory. 
The evolution at small x has no true powerlike behaviour due to the fall of the coupling, but does 
have significant differences from that predicted from a standard NLO in a s treatment. Application 
of the resummed splitting functions with the appropriate coupling constant to an analysis of data, 
i.e. a global fit, is very successful. 
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1. Introduction. 

There has recently been a great deal of interest in the solution to the BFKL equation [1], 
triggered by the calculation of the NLO correction to the kernel [2] [3] and the apparent result that 
this leads to a huge correction to the LO result. A number of subsequent papers have examined the 
solutions to this equation and/or its consequences [4]-[9] drawing a variety of conclusions (dominant 
negative NLO anomalous dimensions, oscillatory behaviour, non-Regge terms, instability, break- 
down of perturbation theory), most being rather pessimistic. This has prompted work on ways to 
at least estimate contributions to the kernel at even higher orders, and obtain perturbative stability 
via a summation [10] [11]. 

I will take the point of view that the most significant result of the NLO kernel is that it 
indicates very strongly how the coupling constant should run in the BFKL equation, i.e. that the 
scale in the coupling should be chosen to be the transverse momentum at the top of the gluon ladder 
k 2 . Making this choice [12] I follow many previous authors in examining how this choice affects the 
solution to the LO equation [13]-[21]. I find that at leading twist the solution factorizes into a part 
dependent of the input to the equation, but independent of the scale k 2 , and a part independent of 
the input which governs the evolution in k 2 [14] [15] [21]. The former is disastrously contaminated by 
the diffusion [22] into the infrared, and without a low k 2 regularization is indeterminate due to the 
presence of infrared renormalons giving behaviour ~ exp(— n/?o(ln(Qo/A 2 )) 3 /(.A 2 ln(l/x))), where 
Qq is the scale of the input to the equation, n is an integer, and A ~ 4. This is entirely consistent 
with Mueller's result [23] [7] on the range of applicability of the BFKL equation. It renders the NLO 
correction to the kernel which is not associated with running of the coupling rather unimportant 
since the infrared contamination renders even the LO result untrustworthy. 

However, the part of the solution governing the evolution in k 2 is not only infrared safe but 
is influenced strongly by diffusion into the ultraviolet. Hence the effective scale in the problem 
is greater than k 2 , and this increase becomes more significant as x decreases. This leads to 
the effective coupling constant decreasing as we go to smaller x, behaving like l/(lnfc 2 /A 2 ) + 
^(ln(l/x)/ln(A; 2 /A 2 ))2) rather than l/(ln(A; 2 /A 2 )). This result is quantified by using the BLM 
scale fixing procedure [24] for both LO and NLO quantities, obtaining precisely the same result 
of A = 3.63 in both cases. It suggests that the effective splitting function governing the evo- 
lution does not grow like a power of x _1_A as x — > 0, but is softened to something of the form 
l 

(1/x) exp((ln(l/x)) 2 )p(/c 2 )), though it seems difficult to obtain the precise form. This result means 
that the NLO corrections to the kernel not concerned with the running of the coupling are also 
relatively unimportant for the term governing the evolution, simply because the coupling constant 
associated with them is so small. Therefore, it seems as though we have good predictive power for 
the evolution of the gluon at small x, but that it is very different from the LO-BFKL prediction 
with fixed a s . Because the behaviour of physical structure functions at small x is related to the 
gluon via the convolution of a fc 2 -dependent cross-section at the top of the gluon ladder [25] [26], all 
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such effects are associated with the ultraviolet diffusion. Hence, the evolution of physical quantities 
is governed by the same effective coupling constant, and is completely predictive, being somewhat 
different from both the LO-BFKL predictions with fixed a s and the fixed order in a s (Q 2 ) DGLAP 
descriptions. 

In this paper I will demonstrate the results discussed above. I will start with a brief discussion 
of the LO BFKL solution with fixed coupling, emphasising the role played by the infrared and 
ultraviolet regions of transverse momentum. I will then look at the same equation for running 
coupling, showing how the solutions change. This will facilitate a discussion of the real importance 
of the total NLO correction to the BFKL equation. Finally, I will examine the implications of my 
results for physical quantities and give a brief outline of phenomenological consequences, showing 
that my results work very well when used to analyse experimental data. I note that a very brief 
account of this work, which nevertheless contains many of the main ideas, appears in [27]. 

2. BFKL Equation for fixed a s . 

The BFKL equation for zero momentum transfer is an integral equation for the 4-point kx~ 
dependent gluon Green's function for forward scattering in the high energy limit, f(ki,k2,a s /N) 
where N is the Mellin conjugate variable to energy. Throughout this paper I will consider the 
canonical physical process of deep-inelastic scattering where the bottom leg is convoluted with a 
bare gluon density and the top leg with an off-shell hard scattering process. Hence, kz is taken to 
be some fixed scale Q\ typical of soft physics 1 , while k\ = k 2 , i.e. a variable scale typically 3> Q^. 
In this case N is the conjugate variable to x, i.e. we define the moment space structure functions 
by the Mellin transformation, 



and the moment space parton distributions as the Mellin transformation of a rescaled parton density, 




(2.1) 



i.e. 




(2.2) 



Using these definitions the BFKL equation becomes 



f(k 2 , Ql a s /N) = f(k 2 ,Q 2 ) + § ^ V, k 2 )f(q 2 ), 

iV Jo 1 



(2.3) 



where 



K°(q 2 ,k 2 )f(q 2 ) = k 2 



f(q 2 )-f(k 2 ) 
| k 2 — q 2 | 



+ 



(V + fc 4 )2 



f(k 2 ) 



i 



) 



(2.4) 



1 Strictly speaking, within the leading twist collinear factorization framework this lower leg should be 
on-shell, so Qq is a regularization scale. 
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f°(k 2 , Qq) is the zeroth order input and a s = (3/tv)cx 8 . As a simple choice I take 

f°(k 2 ,Q 2 ) = 5(k 2 -Q 2 ), (2.5) 

i.e. the incoming gluon has a fixed nonzero virtuality. With this definition a moment space gluon 
structure function can be defined as 2 



Q(Q 2 ,N) = j ^f(N,k 2 ,Q 2 )xg B (N,Q 2 ), (2.6) 



rQ Ah2 
'0 

where gs(N, Qq) is the bare gluon distribution as a function of the factorization scale Qq- 3 

In order to solve this equation it is convenient to take a further Mellin transformation with 
respect to k 2 , i.e. define 

/>oo 

/( 7 ,A0= / dfc 2 (A; 2 )- 1 -'7(fc 2 ,JV)- (2-7) 
J o 

This leads to the BFKL equation written in the form 

7(7, TV) = /°( 7 ,Q 2 ) + (a s /N)x(7)f(j,N), (2.8) 
where /°( 7 , Qq) = exp(— 7 ln(Q 2 )) and x(l) is the characteristic function 

x(7)=W-^7)-W-7). (2-9) 

Hence, 

For asymptotically small x this can be accurately inverted back to x and k 2 space using the 
saddle point technique to give the celebrated result, 

k 2 \* - In 2 (k 2 /Q 2 " 



■-)> (2 - u) 

where A = 41n2a s and ••• denotes subleading terms as x — > 0. Hence, we see that the BFKL 
equation at LO predicts powerlike growth in x~ x and in k 2 , as well as a diffusion in k 2 . One can 
also be a little more systematic and solve for the coefficient functions and anomalous dimensions 



2 In this paper I will ignore the singlet quark distribution. This is purely for simplicity and does not 
change any of the conclusions at all. In most expressions the replacement of gs(N, Qq) with gs{N, Q ) + 
|Sb(JV, Qq), where Eb(A^, Qq) is the bare singlet quark distribution, is all that is required to make them 
completely correct. 

3 In making this definition of the gluon distribution we have defined a factorization scheme. 
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for the gluon, it is easy to generalise (2.10) to give the double Mellin space expression for the gluon 
structure function 

~ f°(>y,N)g B (N,Ql) 

g( 7 ,A/) = — — (2.12) 

7(1 - (a a /N)x{V) 

and 

1 P2+ 100 

Q(Q\N) = — ci 7 exp( 7 ln(Q 2 ))a(7,A r )- (2.13) 

From (2.12) we see that there are poles when 1 — (a s /N)x{j) = 0. Defining the rightmost solution 
of this equation by 

x(VWtf))^ = 1, (2.14) 
we obtain the leading twist solution for the gluon structure function 

"|2 AT\ — ^ ( (AT 



^' -wiww ilJ 3M - (2 - 15) 

Hence, 7 (cL j /A r ) is the anomalous dimension governing the Q 2 evolution of the gluon [28], and 
1Z(a s /N) = — (a s /A r 7 x'(7°)) _1 is a type of coefficient function giving the normalization [25]. Each 
of these may be expanded as power series in (a s /N), which then lead to power series in a s ln(l/x) in 
x-space. Both are only convergent for a s /N < 4 In 2, each developing a branch point showing that 
in x-space they grow like x~W~ A . Using the saddle point technique one may find the asymptotic 
form of the x-space splitting function and coefficient function finding that 

P°{x) - ^*" A ( r i ) , (2-16) 



x 



(56vrC(3))2(a s ln(l/x)) 3 / 2 



and 

1 

-A ' ' > 2 



R(x) = 41n2a s x~ A - - - - - . (2.17) 

V ; Vl4vrC(3)a s ln(l/x); V ' 

Therefore, both the anomalous dimension and the coefficient function predict powerlike behaviour 

for the gluon distribution, although the true input for the distribution is really R(x) convoluted 

with 5b(x,Qq) °f course, and this leads to the exact form of R(x) being sensitive to the input 

f°(N, Qq). 4 However, this powerlike behaviour does not set in until very small x, as may be seen 

by examining the terms in the expansion for each quantity in powers of a s ln(l/x). 

It has long been suspected that the diffusion property of the solution to the BFKL equation 

may have serious consequences when working beyond the strictly LO framework [13] [22] [29] [20] [23]. 



4 In the language of the factorization theorem this translates into R(x) being regularization scheme 
dependent, e.g. if one uses dimensional regularization rather than a off-shell gluon, R(x) has a factor of 



(a s ln(l/x))- 3/4 rather than (a s ln(l/x))~ 2 . 
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One may appreciate this by recognising that in the small x limit, defining £ = ln(l/x), we may 
write, 

i(k 2 ,Qt,0 = J drH{k\r\i')i{r\Qli-i'). (2.18) 

For a given £' we can ask for the mean and the deviation of ln(r 2 ). This is equivalent to asking 
for the typical ln(/c 2 ) at some point along the ladder diagram representing the function f(/c 2 , Qo,£)j 
and also its spread, i.e. the range of important values of q 2 involved in finding the solution of the 
BFKL equation. The result is well known: 



< Hr 2 /(kQ )) >= f 1 - 2 | ) , (2.19) 



and the RMS deviation is 

a 2 



28C(3)a s £'(l-|). (2.20) 



So over much of the ladder < ln(r 2 ) >~ ^ ln(/c 2 /Qq) and a ~ (14£(3)a s ln(l/x))2 . Hence, for 
very low x there will be significant diffusion into both the infrared and the ultraviolet. In the case 
of fixed coupling this does not cause any serious problems. However, in the case of the running 
coupling the size of the coupling grows quickly in the infrared, and hence this diffusion suggests 
that there will be serious contamination from nonperturbative physics. 

Before looking at the BFKL equation for running coupling let us briefly examine the role played 
by the various regions of q 2 in the fixed coupling case. In order to determine the role played by the 
region of low transverse momentum we consider a upper cut-off of k 2 in the integral in (2.3). The 
only restriction we place on k 2 is that kg <C k 2 for whatever transverse momentum we ultimately 
wish to consider at the top of the ladder. With this restriction we see that for all momenta over 
which we integrate we have the relation 

K (q 2 ,k 2 )f(q 2 ) = f(q 2 ) + O^j, (2.21) 

and inserting into the cut-off version of (2.3) we obtain 

^ K°( q 2 ,k 2 )f(q 2 ) = ^h(k 2 ,f(k 2 < fcg)) + o(§). (2-22) 

The integral over the region q 2 < k 2 <C k 2 contributes only a constant to the right hand side of 
(2.3), dependent on the form of f(q 2 , N) at low momentum, but independent of the value of k 2 we 
consider as long as it is large. If k 2 is actually smaller than k 2 then /i(/cq,/(/c 2 < k 2 )) becomes a 
much more sensitive function of k 2 , and in the limit fe 2 — > it is easy to see that it becomes equal 
to the value of the integral in (2.3) with no upper cut-off. Hence, h(k 2 , f(k 2 < k^)) has the same 
structure for k 2 — > as the full integral on the right hand side of (2.3), but tends to a constant 
function of k^ for k 2 » kg. 
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Thus, if we imagine imposing an infrared cut-off on (2.3) we can simply subtract the result of 
the integral up to k 2 , (now with a different f(q 2 ) for low q 2 , in particular the infrared cut-off renders 
it infrared finite) from the right hand side of (2.3), obtaining (up to higher twist corrections). 

f(k 2 , Ql a./N) = f(k 2 ,Q 2 ) - ^h(k 2 , k 2 ) + ^[ ^-K°(q 2 , k 2 )f(q 2 ). (2.23) 

Taking the Mellin transform of this equation we get 

/( 7 ,7V) = /°( 7 ,Qo) " (a./#)M7,*o) + (a s /N) X (j)f(j, N), (2.24) 

where ^(7, fcg) is analytic for 7 > (h(k 2 , k 2 ) tends to a constant at high k 2 ). This second term on 
the right may simply be absorbed into the definition of the input and our expression for 7(7, iV) is 
exactly the same as in (2.10) up to this transformed input, i.e. 

f (7 AO - f°^K)-(a s /N)h( % k 2 ) (2 25) 

Performing the inverse Mellin transformation, then for the leading twist solution the pole is in ex- 
actly the same place and we obtain exactly the same k 2 dependence as previously, but a potentially 
very different iV-dependent normalization. Hence we see that the region of transverse momentum 
<C k 2 contributes very significantly to the overall normalization of our leading twist solution, but 
negligibly to the evolution, essentially because the contribution from the infrared region coming 
from the convolution in the BFKL equation is the same for all high k 2 . We also notice that the 
other, higher twist poles found in (2.10) are now eliminated by the presence of —(a s /N)h(-y,k 2 ). 

This above argument is hardly new, and much more detailed analysis can be found in [30] 
[31] who consider the Mellin space solution carefully, showing that the infrared cut-off does indeed 
change only the residue of the rightmost pole in 7 (and removes all poles in the left half plane). 
It is also noted that infrared cut-offs influence only the normalization of the gluon distribution, 
leaving the shape in x as well as Q 2 largely unchanged [29]. This is because the effect is to change 
the type of singularity in iV-space, but not the actual position, i.e. N = 41n2a s . However, this 
is not usually discussed together with the phenomenon of diffusion. In the case of fixed coupling 
the effect of diffusion is less important than for running coupling for the obvious reason that the 
coupling is the same at all scales. Nevertheless, the above arguments imply that in the case of 
running coupling diffusion into the infrared, i.e. strong coupling, should again only influence the 
normalization of the gluon, while diffusion into the ultraviolet, i.e. weak coupling, should only 
influence the evolution in Q 2 . We will now investigate this in more detail. 
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3. BFKL Equation for Running Coupling. 

It was expected in [12] that the way to incorporate the running coupling into the BFKL 
equation was to modify (2.3) to 

f(k 2 , Ql a s (k 2 )/N) = f°(k 2 , Ql) + MP r ^K°(q 2 , k 2 )f(q 2 ), (3.1) 

iV Jo Q 

where 

as = l/(f3oHk 2 /A 2 )), (3.2) 

/3o = (11 — 2iV//3)/(47r), and Nf is the number of active flavours. One of the main results of the 
NLO corrections to the BFKL kernel is to show that this is indeed an effective way to account 
for the running coupling (this will be discussed more later). One can solve this equation in the 
same type of way as for the fixed coupling case, i.e. take the Mellin transformation with respect to 
(fe 2 /A 2 ). It is most convenient to first multiply through by ln(/c 2 /A 2 ), and then obtain 

where fto = (vr/?o/3). The inclusion of the running coupling has thus completely changed the form 
of our double Mellin space equation, turning it from a simple equality into a first order differential 
equation. However, this may be easily solved to give, 



/( 7) AO = exp(- A( 7 , N)/(0qN)) f°° d ^^Qo) exp(A( 7 )/(A) N))d% (3.4) 



'7 

where 



*(7) = /i 7 x(7)rf7^(2^(l)(7-^)-ln(^^y)). (3.5) 

The leading singularity in the 7 plane for exp(— X{^)/{(3q A)), is cancelled by an integral from 
— > 7 of the integrand depending on 7 [15], and so up to higher twist corrections we may simplify 
(3.4) to 

/(7, A) = exp(- X( 7 )/(A) A)) [°° d /o(7,A,Q 2 ) exp{X (j)/(p N))d^. (3.6) 

Jo «7 

Using our previous choice of input, i.e. fixed virtuality, we obtain the moment space gluon structure 
function 



1 . 

^(Q 2 ,A) = — / -exp(7ln(Q 2 /A 2 )-A(7)/(^oA))d 7 

l-OO 

x / exp(-7ln(Q 2 /A 2 )+A(7)/(^ A))d7to(Qo^) 
J 

= a(Q 2 ,A)&(Q 2 ,A) ffB (Q 2 ,A). 



(3.7) 
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Hence, as in the case of fixed coupling constant, at leading twist the solution has factorized 
into a Q 2 -dependent part a(Q 2 ,N) which determines the evolution, and an input dependent part 
b(Ql,N) which can be combined with the bare input gluon distribution to provide the input for 
the gluon distribution [14] [15] [21]. This time the different parts are not so easy to calculate though. 
Clearly the behaviour of both functions is determined by the form of exp(X(7)/(/3 A r )), since this 
determines the singularity structure. 

Considering b(QQ,N) we find that exp(X (7) / \(3q N)) has poles at all positive integers, and 
zeroes at and all negative integers. Hence, b(Qo,N) is not properly defined, since the integrand 
has an infinite number of poles lying along the line of integration. These are due to the divergence 
of the coupling at low k 2 and can only be removed by some infrared regularization. Hence, the 
diffusion into the infrared has destroyed the apparent (limited) predictive power for the input. 
Imposing some regularization scale k 2 . and repeating the same arguments as the previous section 
it is clear that up to higher twist corrections the effect of the regularization is simply to leave the 
factor a(Q 2 ,N) unchanged, and change b(Ql,N) to 

c(Ql k 2 , N) = (^M^ + h(j, k 2 )^j eMX(l)/(PoN))dj, (3.8) 

where the factor + h{j, k 2 )^j removes the singularities in exp(X( ; y)/(f3 N)). Thus, we 

have 

G(Q 2 ,N) = a(Q 2 ,N)c(Q 2 ,k 2 ,N)g B (Q 2 ,k 2 ,N), (3.9) 

as a well-defined solution. 5 For a given regularization one can solve for c(Qq, k 2 ,N), as has been 
done numerically 6 , generally obtaining some powerlike growth in x-space, but which is totally 
dependent on the type and scale of regularization [16] [17] [20] [32]. No real predictive power remains 
(this will be discussed more in §5). 

Even without regularization there is no obstruction to solving for the Q 2 dependent part of the 
gluon distribution, and this is unchanged by this regularization, i.e. is unaffected, up to 0{k 2 /Q 2 ) 
corrections, by the diffusion into the infrared. The function a(Q 2 ,N) is, of course, determined by 
the singularities of exp(— X (<y) / (flo N)) in the 7 plane. Here we notice a fundamental difference 
between the cases of the fixed and running couplings. Whereas previously the leading singularity 
was a pole at (a s /N)x( , y) = 1, i.e. at 7 — > \ as N — > 41n2a s , now the leading singularity is an 
essential singularity at 7 = 0: there is no powerlike behaviour in Q 2 . Similarly, the branch point in 
the N plane at 4 In 2a s has become an essential singularity at N = 0: there is no powerlike behavior 

5 That the solution at leading twist is of this general form was shown in [21] by putting the BFKL 
equation with running coupling in the form of an infinite order differential equation with effective potential 
depending on the low k 2 regularization of the coupling. 

6 The numerical solutions are always for the whole of the gluon structure function, not just c(Q 2 ,k 2 ,N). 
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in x. The introduction of the running of the coupling has therefore also had an extreme effect 
upon the evolution, changing its character completely. This point has been noticed before [15] [20], 
but not emphasised or studied in detail. Hence I stress the fundamental results of introducing 
a running coupling: the Q 2 -independent part of the solution is formally divergent, and hence is 
totally regularization scheme dependent: the Q 2 -dependent part has no powerlike growth in x. 

In fact we can obtain some information about the x behaviour by noting that we can find the 
inverse Laplace transformation of exp(— X(i)/(J}qN)) precisely [20] [8]. It is a standard result that 

1 r io ° l l 

— / exp(NS + K/N)dN = (A/Z)*I 1 (2(M)S), (3.10) 
2m J_ ioo 

l 

were h(z) is the modified Bessel function, which for large values of its argument — > exp(z)/(2irz) 2 . 
Hence for large £ 

a(£, 7 ) ~ (-X(7)/^ 3/4 exp(2(-^( 7 )/^))i (3.11) 

It is difficult to perform the inverse Mellin transformation to get the Q 2 dependence, but the leading 

singularity is at 7 = 0. Thus, for any Q 2 the leading twist solution for a(£,Q 2 ) must have small 

1 

x behaviour going like exp(£2) rather than the exp(A£) for the fixed coupling case. This is easy 

to understand in terms of the diffusion picture. Since the function a(£, Q 2 ) is insensitive to the 

diffusion towards the infrared, but sensitive to that into the ultraviolet, we expect the typical scale 

in the process to be determined by this latter diffusion. Thus the typical scale for the process will 

1 

be approximately set by ln(A; 2 ) ~ ln(k 2 ) + a ~ ln(k 2 ) + 4(q s (/c 2 )£) 2 . Hence, the effective strength 
of the running coupling will be set by k 2 , rather than k 2 , and as x — > we will have ofjt ~ l/(£)2 . 
This type of effective coupling has precisely the effect of turning the low x behaviour of the fixed 
coupling solution to that which we find for the running coupling. Hence, the diffusion into the 
ultraviolet has a major impact on the Q 2 dependent part of our gluon distribution, but in a well 
controlled, and in principle calculable way, unlike the effect of the infrared diffusion on the Qq 
dependent input. 

Of course, this is just a qualitative argument giving only the general form of the results. It 
is also for the function &(x,Q 2 ), which must be convoluted with an unknown, if (^-independent 
input function in order to obtain physical results. It would be nice to be more quantitative, and 
also to calculate some physical quantity unambiguously. For example, staying in moment space we 
can examine (dQ(N, Q 2 )/dln(Q 2 )), which is an entirely perturbatively calculable quantity, and its 
transformation into x-space. This will be considered in the next section. 
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4. Solving the BFKL Solution for Running Coupling - Evolution. 

The usual approaches taken to finding the solution for a(Q 2 ,N) (or the full solution) are to 
assume that for small x one can expand X(j) about 7 = \ to some finite order in 7, 7 usually to 
C(7 3 ), or to use the saddle point method. Neither of these are at all accurate unless Q 2 is very large 
indeed. This is because along a line parallel to the imaginary axis X(j) is not at all well represented 
by the first few terms in a power series in 7 about either 7 = \ or about the saddle point. The 
former can be seen in fig. 1, where we compare the the full function exp(7 ln(Q 2 /A 2 ) — X(j)/(/3oN)) 
to the case where the exponent is truncated at 0(-y 3 ) along the line 9J( 7 ) = |. 8 Clearly the integral 
over the two contours need bear little similarity. 

When using the saddle point technique one finds the minimum of the exponent of the integrand 
in the definition of a(Q 2 , N) and expands in a Taylor series about this point. This minimum occurs 
when 

^( 7 ln(Q 2 /A 2 ) " X( 7 )/((3 N)) = 0, (4.1) 
which using the definition (3.5) leads to 

X(7) - ^4(7) = 1, (4.2) 



PoHQ 2 /^ 2 )N N 
i.e. at 7 = j°(a s (Q 2 )/N), the anomalous dimension for the fixed coupling case, but with the 
running coupling evaluated at scale Q 2 . The integrand defining a(Q 2 ,N) is thus evaluated along 
the axis K( 7 ) = j°(a s (Q 2 )/N)), i.e. 

a(Q\N) = -Lexp(y Q 7 °(a s (g 2 )/iV)rflng 2 ) 



/' 

j — 



■ exp( 7 ln(Q 2 /A 2 ) + (X( 7o ) - X( 7 ° + j))/((3 N))d r 



7° + 7 

(4.3) 

Letting, 7 — ► — i 7 and expanding about ^y (a s /N) this becomes 

a(Q 2 , N) = i- exp(| Q 7°(« S (9 2 )/A0)^ 2 J~ (Jp + ■ ■ ■) exp( 7 2 x , (7 )/(2^iV) + ■ ■ -)d>y. 

"°° 7 (4-4) 
This is then normally evaluated by ignoring all those parts not explicitly included above, and 
performing the gaussian integral [33] obtaining 9 

a SP (Q 2 ,N) = -exp( C ~f°(a s (q 2 )/N)dln q 2 ). (4.5) 

^(a s (Q 2 )/N)(- X '(Y>(a s (Q^/N)))2 J 



7 This is equivalent to writing the fc 2 -space BFKL equation as an infinite order differential equation and 
truncating at a low order in derivatives, or iterrating the LO solution in the truncated form (2.12) in the 
NLO equation. 

8 Actually I plot the real part of the functions. The imaginary part is odd and integrates to zero. 

9 A factor of {(3 N/2tt)^ is absorbed into b(Ql, N). 
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This is of the same form as (2.15), i.e. an evolution term governed by the previous anomalous 
dimension and a coefficient function which is a power series in a s /N, except that now a s runs with 
Q 2 rather than being fixed. This could be taken to imply that one can simply extract anomalous 
dimensions and coefficient functions from this solution and that the appropriate scale to use for 
the coupling is Q 2 . 

The invalidity of this assumption is related to the fact that (4.5) is in fact a very poor ap- 
proximation to the full solution for a(Q 2 ,N). This is clear because in x-space both the perceived 
anomalous dimension and coefficient function above grow like xC- 1 )- 411120 ^© 2 ) as x _ > o, whereas 
we know that the complete solution for a(Q 2 ,x) has no real powerlike behaviour in x. We can 
see how we have obtained such a poor approximation by using the saddle point technique if we 
examine the form of the complete integrand along our contour of integration compared with the 
function we have actually integrated making the approximation in (4.4). This is seen in fig. 2, 10 
and it is glaringly obvious that the saddle point estimate is not at all reliable in this case. Formally 
the corrections ignored in evaluating (4.4) are of higher order in (3oa s (Q 2 ) than the terms calcu- 
lated, but their coefficients grow quickly, i.e. like factorials, and to be precise they are powers of 
Po(Xs(Q 2 )(&s(Q 2 )O r higher than the presented results, where r is a positive number, and are thus 
dominant for low enough x. Hence a resummation is really necessary for a true understanding. 

However, an alternative view of the result in (4.5) may lead us towards the correct physics. 
It is not really useful to interpret the prefactor in this equation as a coefficient function which 
tells us something about the normalization of the gluon structure function since a(N, Q 2 ) must be 
multiplied by an unknown iV-dependent function in order to obtain this distribution. Rather, it is 
better to acknowledge that the only real information contained in a(N,Q 2 ) is on the evolution of 
the structure function, i.e. 



dln^(iV,Q 2 ) din a{N,Q 



2-- 



d\n(Q 2 ) dln(Q 2 



= T(N,Q 2 ). (4.6) 



Thus, using a(N,Q 2 ) in (4.6) gives us an entirely perturbative effective anomalous dimension 
governing the evolution of the gluon distribution. Using (4.5) we obtain 

T(A r , Q 2 ) = J° (a s (Q 2 ) /N) - (3 a s (Q 2 ) ^^]~^~~y +^((^ oa *W 2 )) 2 ) r (" s (^ 2 )/^)- 

(4.7) 

So within the framework of the LO BFKL equation with running coupling our unambiguous effective 
anomalous dimension is the naive leading order result with coupling at scale Q 2 plus a series of 
corrections going like powers of (3oa s (Q 2 ). 

It is tempting to interpret the whole solution for T(a s (Q 2 ), N) as simply telling us the appro- 
priate scale to use in the coupling constant for the normal LO result. Indeed, this is the philosophy 
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Again I plot only the real part. 
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in the BLM scheme [24] for scale fixing which uses the NLO /^-dependent corrections for any pro- 
cess to determine the scale to use for the coupling in the LO expression. However, in this case of 
an anomalous dimension for a structure function we have to decide whether it is appropriate to do 
this in iV-space or x-space, i.e. should we write 

* r(N,a s (s(N)Q 2 ))G(N,Q 2 ), (4.8) 

or i 

^Inm * i P (z,aMz),Q 2 ))G(x/z,Q 2 )dz. (4.9) 

Since the moment space expressions are less physical, being defined only by analytic continuations 
over much of the iV-plane we choose the latter. 11 As we will see later, this decision is backed up 
by higher order calculations. Note that both T(N, Q 2 ) and P(x,Q 2 ) are entirely independent of 
factorization scale, and are functions only of renormalization scale. Indeed, if there were a direct 
probe of the gluon, i.e. G(x, Q 2 ) were directly measurable, then both T(N, Q 2 ) and P(x, Q 2 ) would 
be physically defined quantities. As such the choice of the renormalization scale is entirely open. 

The simplest thing we can do is to choose the scale for the coupling constant in the leading 
order expression so that the NLO term in the x-space version of (4.7) is exactly produced by the 
expansion about a s (Q 2 ). Thus, writing this x-space expression as 

(x/a s (Q 2 ))P(x, Q 2 ) = P °(a s (Q 2 )0 - o a s (Q 2 )p\a s (Q 2 )O + O((P a s (Q 2 )) 2 )r(a s (Q 2 )Z), (4.10) 

this is the same as 

(x/a s (Q 2 s(Ca s (Q 2 ))))P(x, Q 2 ) = p°(a s {Q 2 s(^a s {Q 2 ))) + O((/3 a s (Q 2 )) 2 )r(« s (Q 2 )e), (4.11) 
if we choose 

This is the usual BLM scale fixing, but here we have extra information since, in principle at least, 
we know higher order terms and we would expect r{a s {Q 2 )^) to be small if the scale fixing is 
correct. 

Equation (4.12) can be solved for arbitrary x, but it is first useful to examine the limit of 
x — > in order to see if our previous expectations based on qualitative arguments are confirmed. 
Hence we need each of the terms in (4.12) in this limit. As x — > 

p°(a s (Q 2 )0 - r exp(A(Q 2 )0(a s (g 2 )e)" 3/2 , (4.13) 

(56ttC(3))2 



11 Fixing the scale in iV-space would lead to a scale which was singular at N = A(Q 2 ), which does not 
seem a sensible proposition, while in s-space it is a smooth function of x as we will see. 
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and therefore 

d\n(a (Q*)) ~* rexp(AO(« s (Q 2 )0 2. (4.14) 

In order to calculate the x — > limit of j5 1 (q s (Q 2 )0 it is easiest to first consider its moment space 

analogue, i.e. the second term on the right of (4.7). First we note that using (2.14) 

d^(a s {Q 2 )/N) x(l\as{Q 2 )/N)) 

dlnMQ*)) X'(7°(a s (Q 2 )M0)' 1 ' ' 

In the limit x — > 0, x(7°) ~~ * 4 In 2 and 7 — > |, but to be more precise, 

x(7 0) ^ 41n2 - 14C(3)(| - 7 °) 2 + ■ • • • (4.16) 

Therefore, 

X'(7°) - -28C(3)(| - 7°) = -28C(3)5 7 °. (4.17) 

Hence, 

d 7 ° X (7°) ^ 2 



dln(a 8 ) X'(7°) ' 7C(3),57 ' (418) 
Since £7° is vanishingly small as x — > we see that the I/7 term in (4.7) becomes subleading 
to the x(7°)/Y(7 ) term. The iV-space version of ^(a^Q 2 )^) is thus (x(7 )x"(7 )/2(x'(7 ) 2 ). 
To progress further we need 5j° as a function of N. This can be obtained by solving (2.14) using 
(4.16). This gives 

This can be substituted into the moment space analogue of p 1 (a s (Q 2 )^) and the inverse transfor- 
mation performed to give 

p\a s (Q 2 )0 - ln2exp(A(Q 2 )0. (4.20) 

This now makes it trivial to solve (4.12), and we find that in the coupling in our LO splitting 
function 

1 

ln(Q 2 /A 2 ) - ln(Q 2 /A 2 ) = ln(Q 2 /A 2 ) + ( 56 ^)?0 2 (a s (Q 2 )^ • (4.21) 

This is exactly the sort of scale change we would expect from the diffusion into the ultraviolet. It 

1 

also leads to xP(a s (Q 2 ), x) ~ exp(1.14(£/a s (Q 2 )) 2 ) as x — > 0, precisely the sort of behaviour we 
would expect from the qualitative discussions in the last section. 

We can also solve (4.12) exactly rather than relying on asymptotic limits using the power series 
expansions of ^ ln ^ s ^2)) anc ^ P 1 (&s{Q 2 )Q in a s (Q 2 )^. The results of such solutions are shown in 
fig. 3, where I plot the effective coupling constant for Nf = 4 derived as a function of x compared 
to its constant value taking Q 2 as the scale. The qualitative result is entirely consistent with (4.21) 
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though the effective scale is a little smaller than this asymptotic result at highish x values due to 
p°(a s (Q 2 )0 and p 1 (a s (Q 2 )0 not yet having reached their asymptotic expressions. 

Hence, this BLM scale fixing procedure leads to a choice of scale which fits in well with our 
naive expectations, and must be at least broadly qualitatively correct since it does destroy the 
powerlike behaviour we get from fixed order calculations in a s (Q 2 ). Ignoring for the moment the 
fact that we have assumed the manner in which to take account of running coupling effects in the 
BFKL equation (we will discuss possible corrections later), we would still like to know whether 
our prescription is a true representation of the full effect of the running coupling, i.e. whether 
r(a s (Q 2 )£) in (4.11) is really small. At each order in /3oa s (Q 2 ) it is possible to calculate the 
leading behaviour in the limit x — > 0. By power counting one can see that these leading terms come 
from keeping only the next term not explicitly shown in the exponential in (4.4), i.e. the leading 
behaviour is given by 

a(Q\N) = ^expjy 7 °(a s (g 2 )/iV))dW J°° exp^)"^ 7 V VH^hV' V)))^- 

(4.22) 

Expanding the exponential in and performing each relevant integral gives the most singular 
behaviour in (N — X(Q 2 )) at each order in (3 a s (Q 2 ). This is a series of the form 

oo 

a SP (Q 2 ,N) J2(-l) n L n (Poa s (Q 2 )n\(Q 2 )/(N - X(Q 2 ))) 3n/2 , (4.23) 

n=0 

where asymptotically L n ~ (1.84) n n!. Inserting this into (4.6) and performing the transformation 
into x-space leads to a power series of the form 

oo i 

p (x,a s (Q 2 ))^(-l)^ n (/? a s (Q 2 )r(A(Q 2 )a|(Q 2 K 3/2 ) n , (4.24) 

n=0 

where the A n are all positive. If A n were equal to 3.63 n /n! the above series would simply be 
exp(— (A(Q 2 )3.63/?oa^ 2 (Q 2 )£ 3 ^ 2 )> which would be precisely the leading correction in the exponent 
of p°(x, a s (Q 2 )) expected using my choice of scale, i.e. 

exp(A(Q 2 )£) -> exp(A(Q 2 )£ - A(Q 2 )e(3.63/3 « 3/2 (Q 2 )^) + •••)■ (4-25) 

In practice this works reasonably well. A\ = 3.63 of course, since this set our scale. A 2 = 7.08 

rather than 6.59, and the terms then slowly increase above (3.63)™ /n\. As n — > oo, A n+1 /A n — > 

1.67/n2, and therefore (4.24) cannot be precisely of the suggested form. Nevertheless, it defines a 

l 

convergent series in (l3oa s (Q 2 )) n (X(Q 2 )a 2 (Q 2 )£ 3 / 2 ) which for a wide range of values mimics the 
desired exponential exp(-\(Q 2 )£(3.63(3 a 3 J 2 {Q 2 )^ ) well. 

This above check is not really terribly useful since the right hand side of (4.25) hardly matches 
exp(A(Q 2 )£) well for very large £, and many other terms are important at all £. Including our scale 
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choice in the LO expression for the splitting function also leads to terms not explicitly shown in 
(4.25) (and in the expansion of the unexponentiated terms in p°(x, a s (Q 2 ))) which are subleading 
in £ at each power of Poa s (Q 2 ) to those discussed above. There are also terms of this type generated 
by the subasymptotic corrections to (4.21). In principle one could compare with terms generated 
form a more careful solution of (4.4), including also the nonleading parts coming from (4.22). This 

rapidly becomes extremely complicated indeed. It appears as though the logarithm of the splitting 

l 

function is indeed an oscillating power series in 

Poa s (Q 2 )(a s (Q 2 )0 2 , but ^ is difficult to prove 
this rigorously. (We do know that the series will converge, or at least be unambiguously summable, 
since the integral defining a(Q 2 ,N) is well defined.) The best check to be done at the moment is 
to calculate the whole of the 0(((3oa s (Q 2 )) 2 ) contribution to the splitting function exactly, and 
compare this to that expected if the scale choice is correct, i.e. 

i (R n <r*x? ( d 2 P°(x,a s (Q 2 )) dp°(x,a s (Q 2 )) \ f p\x,a s (Q 2 )) \ 2 

2 KPo ,W )) y d{lnasm) 2 + 9 (lna s (Q 2 )) ) \dp»{x,a s {Q 2 ))/d(\na s {Q 2 )) ) ' 1 ' ) 

The relevant terms in a(N,Q 2 ) can be found by considering the terms in (4.4) multiplying the 
gaussian which go like 7 6 /{(3qN) 2 and 7 4 /(/3oiV), performing the gaussian integrals and using the 
equality N = a s (Q 2 )j°. This gives 

(n 2 m SP (n 2 m f, R ftf J WBj) r X"(7°) X'"(l°) s X(7°) \ , 

a(Q ,N) = a (Q , N)yl- (3 a s (Q )^_^^-(_ 

(4.27) 

Inserting into (4.6) and making the transformation to x-space we obtain the required 
0{(f3oa s {Q 2 )) 2 ) splitting function p 2 (a s (Q 2 )^). This is compared to (4.26) in fig. 4, where each 
term is weighted by (exp(A(Q 2 )£)) _1 , and the upper limit of a s (Q 2 )^ = 4 is chosen since the first 
20 terms in the series expansions of each expression give a very accurate representation up to this 
value and it easily covers the range relevant for comparison to HERA data. As one can see, above 
«s(Q 2 )£ = 1 the ansatz for the 0(((3oa s (Q 2 )) 2 ) contribution of the splitting function matches ex- 
tremely well to the explicitly calculated value. Below a s (Q 2 )^ = 1 the matching is not so good, 
but this is relatively unimportant since in this region this contribution to the total splitting func- 
tion is small compared to the more leading contributions, i.e. the scale change is quite small and 
P°(a s (Q 2 )£) » f3 a s (Q 2 )p 1 (a s (Q 2 )O > (f3oa s (Q 2 )) 2 p 2 (a s (Q 2 )£). In this region the scale choice is 
also sensitive to the interference with the finite x effects at fixed order in a s (Q 2 ) which are ignored 
using this expansion scheme. Hence, it seems reasonable to conclude that explicit checks strongly 
support the assumption that all the running coupling effects in the evolution can be accurately 
described by the use of the effective scale obtained by solving (4.12) in the LO effective splitting 
function. 
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5. Solving the BFKL Equation for Running Coupling - Input. 

We could also attempt to evaluate b(QQ,N) in the same manner, i.e. expanding about 
7 (<5 s (Qo)/iV) an d calculating an order by order series in /SoQ^Qo)- Of course, without an in- 
frared regulator we know that b{Qfo,N) must be divergent because the integrand has singularities 
along the contour of integration, i.e. at integer values of 7, which lead to ambiguities of order 
(A 2 /(5o) ra , i.e. higher twist. These singularities do indeed show up in this power series solution. 
Expanding about 7 (a s (C t ) 2 1 )/./V) and only keeping the lowest order terms one obtains a sensible 
solution, i.e. 

b SP (Q 2 ,N) = -expf / Q \°(a s (q 2 )/N)dlnq 2 ). (5.1) 

(a s (Q 2 )/N)(-x'(^(a s (Ql)/N)))2 V ) 

Going beyond this approximation one obtains the same sort of series as for a(N,Q 2 ), except that 
because the contour is now along the real axis, rather than parallel to the imaginary axis, the terms 
in the series are all of the same sign rather than oscillating. This leads to at least one power series 
behaviour of the form 

00 

B n n\{faa s {Q 2 )Y(a s (Q 2 )/(N - \{Q 2 ))Y' 2 , (5.2) 

n = Q 

where B n is roughly B n , and B ~ 3.6. One can take the inverse transformation of this series term 

by term, obtaining a power series in x-space which sums to approximately the form 

00 1 
eMHQDO £(BA,a*/ 2 (Qg)£2)». (5.3) 

n=0 

Hence, in this case the power series is suggestive of the fact that due to the diffusion into the 
infrared the appropriate coupling for b(Ql,N) depends not simply on ln((3o/A 2 ) but on ln(Qo/A 2 ) — 
3.63/?o(o!s(<3o)£) 2 i ^e exact opposite of the case for a(Q 2 ,N). 

Of course, the infrared diffusion is a rather complicated problem, and the series in (5.3) is 
only convergent for (B/3 a S s ^ 2 (Qq)£2 ) < 1. This indicates that I have been too simplistic in 
transforming (5.2) to x-space term by term. The series in (5.2) is not defined for any N, and 
before going to x-space we must solve this problem. The series (5.2) may be summed using 
standard Borel transformation techniques. This leads to a well-defined series up to an ambi- 
guity of the form exp(— (N — A(Qq)) 2 " / (B Po&V 2 (Qq)) ■ Now performing the transformation to 

0/9 o — 

x-space we obtain a well-behaved series in 

(B/W (Q 2 ,)^) as well as an ambiguity of order 
exp(A(<3o)£) ex P(—l/(B 2 Po a l(Qo)Q)i wn ere B 2 « 13. This latter ambiguity is due to the pres- 
ence of an infrared renormalon [34] in the expression for b(QQ, N), and will be cancelled by similar 
ambiguities in higher twist corrections. 12 Such terms are therefore taken to estimate the size of 

12 The ambiguity is seen as the nonperturbative contribution to the solution in [9]. 
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higher twist effects. In this case we see that due to diffusion becoming enhanced at small x, this 
infrared induced uncertainty quickly becomes large at small x, and indeed the calculation of the 
normalization of the gluon Green's function is only at all reliable in the limit 

13f3 2 a 3 s (QDt « 1. (5.4) 

Hence, we find that even if we had a reliable model for the bare gluon distribution gsiQl, N) 13 we 
cannot calculate the input for the gluon distribution at small x within perturbation theory, and 
previous conclusions on the infrared diffusion physics ruining perturbative predictability [23] [7] are 
confirmed. In particular we note that the requirement in (5.4) is basically identical to that found 
in [23] [7], and indeed, if the series in eq. (45) of the latter is summed it has an ambiguity of exactly 
the same type as discussed above (though in [7] the series in x-space was found directly). However, 
here I stress that this ambiguity is unique to the normalization function, and does not affect the 
evolution, which is calculable in perturbation theory. 

Before finishing this section it is interesting to discuss the relationship between the solutions 
obtained via the techniques in this paper and solutions obtained by other authors. It has been noted 
by several authors ([7] [8] [9]) that the asymptotic solution for the BFKL equation with running 
coupling has the general form, 

f(Q 2 , Ql ~ exp (\(QQ )£ + KPlaie ~ +-), (5-5) 

V 56C(3)a s £ ; 

where unless explicitly stated a s is at some fixed scale fi, and K = (7/6)£(3)(3/7r) 3 (41n2) 2 . This 
seems rather at odds with the results discussed above. However, it is not difficult to see how this 
solution appears. Ignoring the term in the exponent going like a^£ 3 one achieves a solution of this 
form simply by taking the transformation to x-space of the product a SP (Q 2 , N)b SP (Qq, N) in the 
limit x — > 0, and only keeping the most dominant terms in the series expansions of the couplings 
about scale fx. 

It is not too much extra work to see where the a^ 3 terms come from. Consider if rather 
than taking the saddle point approximations for a(Q 2 ,N) and b(QQ,N) one takes the solution of 
(4.22) for a(Q 2 ,N) and the equivalent expansion for b(Ql,N). The solution for b(Ql,N) in this 
approximation is of precisely the same form as (4.23) once we replace a SP (Q 2 , N) with b SP (Ql,N) 
and remove the factors of (—1)™ (the L n are identical). If we multiply the two series in these 
expressions together then since at large orders L n ~ a n nl, the resulting series is to good accuracy 
proportional to 



oo 



EW(/W(A/(iV-A)) 3 r, (5.6) 

n=0 



13 Given that the function b(Qo,N) is dependent on the type of collinear regularization as well as the 
ambiguity discussed above this actually seems rather unlikely. 
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where a s = a s (fi 2 ), i.e. we expand cx s (Q 2 (Qq)) about a s (fi 2 ), and asymptotically L 2n +2/^2ra — * 
(63C(3)/(81n 2))n 2 . Multiplying this by the two saddle point solutions, and performing the trans- 
formation to x-space this sum introduces precisely exp(K j3 2 a^ 3 ) with the correct value of K. 
Hence, this non-Regge term comes about due to interference between the input term b(Ql,N) and 
the evolution term a(Q 2 , N). 

Hence, these previous results do appear by taking the transformation to x-space of the product 
of truncated solutions for a(Q 2 , N) and b(Ql, N). However, I would argue that these solutions are 
not representative of any real physics, since neither of these truncations is at all accurate except at 
quite high x. For given Ql (5.5) is only applicable for x satisfying (5.4), in which case the x —> 
approximations used to derive (5.5) are generally rather inaccurate. When (5.4) is not satisfied 
the transformation of b(Ql,N) is indeterminate, and that of a(Q 2 ,N) requires resummation. The 
only sensible option seems to be to factor out b(Ql,N) and simply use a(Q 2 ,N) to determine the 
evolution as accurately as possible, rather than trying to find f(Q 2 , Qq, £). Then we know from the 
general arguments already discussed that the Regge term exp(A£) is nothing to do with the true 
result, let alone the non-Regge term ex.p(K/3 2 a^ 3 ). 

6. NLO Corrections to the BFKL Equation. 

So far I have simply assumed that an accurate way to account for the running of the coupling 
in the LO BFKL equation is to use (3.1). This is an assumption which involves the resummation of 
an infinite number of terms, i.e. it assumes that at all orders in a s (/j, 2 ) the dominant contribution 
to the BFKL equation due to the running coupling is 

§(-ir(/W^ 2 ) Hk 2 /f)Y r ^f*V, k 2 )f(q 2 ), (6.1) 
JV Jo Q 

Until recently this has been an assumption for all n > 1 although the above terms must be present. 
However, the recent calculation of the NLO correction to the BFKL equation has given us some 
insight into this question. Formally the NLO BFKL equation may be written as 

f(k 2 ,Q 2 ,a s ^ 2 )/N) =f(k 2 ,Q 2 ) + ffip-} f^-(K°(q 2 ,k 2 ) 

- f3 a s {^ 2 ) Hk 2 /^)K°(q 2 , k 2 ) - a s (p 2 )K\q 2 , k 2 ))f(q 2 ), 

(6.2) 

where K 1 (q 2 , k 2 ) can be found in [2]. This is the strictly NLO equation with no resummation at 
all. The separation of the NLO part into the running coupling part and the part depending on 
K 1 (q 2 , k 2 ) is arbitrary. The former is the first term in the infinite series we have already considered, 
but the latter also contains some pieces which may be associated with the running of the coupling, 
i.e. going like (3q. 
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This equation can be solved using the same methods which were applied in §4. Taking the 
Mellin transformation, this time with respect to (k 2 / ft 2 ) we obtain 

/(7, N) = /°( 7 , Ql) + {^p-) ((%o(7) " « s (^ 2 )Xi(7))/(7, N) + (3 a^ 2 ) d ^o(7)/(7)) ^ > (6 3) 

where x 1 ^) can & l so be found in [2]. As in §3, this is a first order differential equation in 7, and it 
can be solved in the same manner. In fact it is rather easier to alter (6.2) to the slightly different 
form 

f(k 2 ,Q 2 ,^ 2 ) = f(k 2 ,Q 2 ) + jf ^-(K\q 2 ,k 2 ) - a s ^ 2 )K\ q 2 ,k 2 ))f( q 2 ). (6.4) 

This is identical to (6.2) up to NNLO in a s (fi 2 ) and is a common way for the NLO BFKL equation 
to be written since it makes the solution easier. One must simply remember that the solution 
obtained is only uniquely defined up to NLO in a s (fi 2 ) when the coupling a s (k 2 ) is expanded 
about a s (fi 2 ). If we take the Mellin transformation of (6.4) with respect to (k 2 /A 2 ) we obtain 

df(>y,N) _ d/°( 7 ,Qg) 1 , , , 



dj dj (5 N 



(Xo(7)-« s (^)Xi(7))/(7,A r ), (6.5) 



which is identical to (3.3) except for the NLO in a s (p> 2 ) correction to the kernel. It can therefore 
be solved in exactly the same manner as this previous equation (this would also be true for (6.3)), 
again obtaining a solution factorizing into a <5 2 -dependent part and a Q$ dependent part. Each of 
these is a contour integral and analogously to the previous treatment expanding about the saddle 
point when performing the inverse Mellin transformation to Q 2 - or Qg-space produces an ordered 
series in a s (/i 2 ), as long as we also expand 0! s (Q 2 (Qq)) about a s (/j, 2 ). This time the saddle point 
is at [35] 

y""'<*./W> = VWW) - A* MQHQIVS-)^^ - a. "ffl^Z +•". <«•«> 

<9(ln(a s )) -Xo(T(as/N)) 

where 

as -x> m* s /N)) = a ^ (Qs/iV) ' (6 - 7) 

is often called the NLO-BFKL anomalous dimension, and all other corrections are beyond NLO in 

Using the previous choice of input we can evaluate the two inverse transformations about the 
saddle point (we only need go further than the strict saddle point approximation when considering 
the (I/7) factor in the Q 2 -dependent integrand to obtain all results up to NLO accuracy - i.e. we 
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For the remainder of this section unless the argument is explicitly stated a s = a s (fi 2 ). 
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use 1/(7° + 7) 1 = I/7 — 7/(7° ) 2 )- This gives a solution for the gluon structure function of the 
form 

G{Q\N) 1 / / 7 1 , ,1 Xo'(7°) , X'i(7°) a x(, / (7°)Xo(7° 



^(Q 2 ,,^) («s/iV)7°(-x[,(7 )) V S V 7° 7 -4(7°) -X[,(7 ) P ° 2 7 °(-4(7 )) 2 
-/Wn(Q //, )Mn(^i-^ + 2 ^C?)J "^-Wo/M )( v - 1 + 2 g ln(as) - x /( 7 o) jj 

x exp / (7 -/^a^lnO? 2 /^ 2 ) — — - — rr- - a^ 1 ) d In q 2 . 
\JQI d(ln(a s )) J 

(6.8) 

This allows us to determine the gluon coefficient function and gluon anomalous dimension up to 
NLO in a s (/j, 2 ), where the former may be defined as the value of (6.8) when Qq = Q 2 , and the 
latter is then determined by the evolution of (6.8) with respect to Q 2 once the coefficient function 
has been subtracted out, i.e. 

-nNLO(- /Ar/o2/,,2\ \ (-> n ( 7"*" , „,1 XoCt) , Xl (l° 

U (a s /N,Q In ) = l- a J- + 7 + 



(«s/iV)7 (-X^,(7°))V V 7° ' -X'oil ) -X'oii 

I *g(7°)Xo(7°) ln(0 2 /a 2 )( 1 I ^ ( 1 + ^ 



27°(-X^(7°)) 2 ^ V 51n(a s )V 7° -x{,(7° 



(6.9) 



and 

7' 



W",«W) ^'-*»-'« ! ^) 5 p|i r «.T' + (-l + l^||l5). (6-10) 



It is gratifying, though necessary, that in both cases the LO results from (2.15) are are repro- 
duced, and the terms ~ (5oa s ln(Q 2 / /u, 2 ) are consistent with the renormalization group. (Note that 
— a s ^ 1 (a s /N) is not really the NLO correction to the anomalous dimension in this scheme - it is 
actually quite similar to the MS factorization scheme anomalous dimension.) Nevertheless, both of 
these quantities are dependent on our choice of input and factorization scheme, and do not contain 
any real physics. 

The only physically unambiguous quantity which may be extracted is the effective anomalous 
dimension defined by (4.6): 

The second term on the right corresponds to the NLO in a s contributions previously accounted for 
when considering the running coupling, while the third gives the additional NLO corrections. By 
examining the part of 7 1 which depends on /3 we can check whether at NLO at least the previous 
assumption about the manner in which to treat running coupling effects was correct, i.e. can see 
whether these do give the dominant contribution at NLO, or whether the conformal parts of 7 1 are 
more important. 
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One can study the terms in (6.11) by finding the explicit form of each as a power series in 
a s /N. However, in the small x limit we can examine the form of the singularities in the iV-plane, 
i.e. the limit of each of the terms as 7 — > \ and N — > A. Using the well publicized fact that 
Xi(^) = 4 In 2 x 6.3 for 4 flavors and in the MS renormalization scheme, and taking the inverse 
transformation back to x-space of (6.11), we obtain 

xP(x, Q 2 ) = a s exp(AO (7^^ - A)« s ( ( ln(Q 2 /^ 2 ) + O.69) - a s f^V) ) ■ (6-12) 

Hence, the last term, although numerically large, is subleading to the effects due to the running of 

1 

the coupling we have previously considered, being a power of (<5 S £) 2 smaller. However, now we can 
be a little more systematic. Examining the full NLO correction Xi(j)> presented in eq. (14) in [2], 
we see that there are contributions which may be interpreted as being due to the running of the 
coupling. These are ^/?o(x 2 (7) + x'(t)) an d — (5/3)/?ox(7)> coming from the NLO correction to the 
reggeon-reggeon-gluon vertex and the purely virtual terms respectively. 15 We imagine that these 
should be moved out of 7 1 in (6.11) and put into the /^-dependent part of the NLO correction. 
Doing this changes (6.12) to 

„, 0.068 a (( 0.188 ln(Q 2 /n 2 ) -0.05 \ n \ ( 1.23 \\ 
xP(x, Q 2 ) = a s exp(A£) - (3 a s V ^ 1 V + °- 69 " <M r ■ 

(6.13) 

Therefore, not only is this additional NLO correction due to the running of the coupling numerically 
very small, but it is also subleading at small x to the terms we have already considered. 16 Choosing 
the renormalization scale [i by setting the (3$ dependent term to zero 1 , we obtain a very minor 
correction to our previous choice of scale for the limit x — > 0, i.e. 

1 

ln(Q 2 /A 2 ) - ln(Q 2 /A 2 ) - 0.26 + 2 (a s (Q 2 )^ , (6.14) 

where in fact there should really an additional constant on the right in the above equation due to 
subleading corrections in (4.21) that I have ignored. The constant on the right of (6.14) is also 
renormalization scheme dependent, though the dominant 3.Q3(a s (Q 2 )^) 2 term is not. 

We can also solve the equation for the scale exactly rather than in the small x limit. Putting 
our additional terms into the definition of the running coupling dependent NLO splitting function, 



15 It does not seem certain whether or not the second of these terms should be included as a running 
coupling effect or not. As will become clear below this is only relevent for the scale choice at high x where 
other considerations from large x terms come into play also. 

16 Not including the — (5 / 3) (3ox(l) term would simply lead to —0.05 becoming 0.26 and 1.23 subsequently 
becoming 0.92. 

17 I choose a s to be a s (Q 2 ) rather than a s (/U 2 ) when doing this. The two are of course equivalent up to 
higher order corrections, but the results of previous sections suggest that this is the appropriate choice. 
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we compare with the previous p l {a s {Q 2 )£) in fig. 5. We see that indeed the corrected p 1 {a s {Q 2 )£) 
is slightly smaller than the original for a s (Q 2 )^ > 1 but is different at higher x, implying a different 
scale choice here to that in §4. Of course, at these higher values of x the differences are not 
too important since, as already mentioned, the scale changes are small here, and there will be 
interference with other effects from the order by order in a s expansion. 

Hence, we find that at NLO our previous assumption about the — Poa s (/j, 2 ) ln(fe 2 /// 2 ) term 
(which had to be present) being the dominant contribution associated with the running coupling 
is very well justified. This gives us confidence, if not a proof, that the approach taken in the 
previous sections, i.e. that the (— /3oa s (/" 2 ) ki(k 2 / [i 2 )) n terms are the dominant contribution from 
the running of the coupling at all orders is roughly correct. Consequently, this full NLO result also 
supports the hypothesis that the LO running coupling effects can be taken account of simply by 
using the x-dependent scale choice, determined by the BLM prescription, in the LO expression for 
the effective splitting function. 

Before considering the details of the NLO corrections to the kernel which are not associated 
with the running of the coupling let us reconsider the NLO BFKL equation. Given the above 
results it seems very unlikely that the NLO BFKL equation as written in (6.4) is will be a good 
representation of the real physics since the overall power of the coupling is allowed to run with 
k 2 while that associated with the NLO kernel is fixed at /j, 2 . Bearing in mind that letting the 
coupling run in the LO equation leads to such dramatic effects, and that at higher orders there 
will definitely be the logs in (k 2 / fi 2 ) associated with the running of this additional factor of a s (/j, 2 ) 
(with what now seem likely to be small corrections) it seems most appropriate to write the NLO 
BFKL equation with running coupling as 

f(k\ Qlf) = f(k 2 ,Q 2 ) + ) jH ^(K°(q 2 , k 2 ) - a s (k 2 )K\q 2 , k 2 ))f(q 2 ), (6.15) 

if attempting to find a complete solution, as proposed in [11]. Strictly speaking a s (k 2 ) should then 
be the two-loop running coupling, but this will make the equation very complicated. I will just use 
the one-loop coupling which leads to a 2nd order differential equation in 7-space 

d 7 2 d 7 2 P N d 1 3/3 2 N Xl{1)J[1, 1 1 } 

This can be solved in a very similar way to the approach in §4, i.e. at leading twist it factorizes 
into the same form as (3.7) 

G NLO (N,Q 2 ) =a NLO (Q 2 ,N)b NLO (Q 2 ,N)g B (Q 2 ,N), (6.17) 

where 

1,. 

1 /"2+* 00 1 

a NLO (Q 2 ,N) = — -exp( 7 ln(Q 2 /A 2 ) - X NLO (j, N)/(p N))d 7 . (6.18) 
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However, X^^od, N) is rather more complicated than the previous X^) It can be expressed in 
the form 

X NLO (l,N)= r X NLo(l,N)di, (6.19) 

where XNLoil-, N) can be written as a power series in N beginning at zeroth order with Xo(t)- As 
seen in [11], though here ignoring any resummations in N, the explicit form is 

Wlo (7, AO = Xo(7) - N^- + * (- (*<2>) 2 - ft (2>W) )+ ... . (0.20) 

Xo(7) Xo V \Xo(7)/ \Xo(7)/ / 

where ^2(7) would also appear at order iV 2 if I had included it. 

It is now possible to obtain some general and rather specific results using (6.16). Putting (6.20) 
into (6.18) we note that the leading singularities in 7 and N are still both at 0, and thus there is 
still no true powerlike growth. Furthermore, the singularity at N = is not affected by any of the 
additional terms in (6.20) beyond Xo(l) since in the exponent in (6.18) the O(N) term leads to 
a constant as N — > and all higher order terms vanish in this limit. Hence, none of these terms 
should affect the solutions in the limit x —> 0, except that the O(N) term should affect the overall 
normalization, and we still expect small x solutions ~ exp((£)2) with the exponent the same as 
in the LO case. Hence, higher order corrections to the BFKL equation should be very subleading 
when calculating physical quantities. This implies that the scale for the coupling in higher order 
corrections should be of the same type as at LO, i.e. falling with x. 

It is also possible to be more quantitative. (6.18) can be solved using the same techniques as 
in §4 - expanding about the saddle point leads to an ordered expansion in a s (Q 2 ). Using (6.20) it 
is easy to find that the saddle point is now at 

j SPNLO (a s (Qi)/N) = 7 °(« s (Q 2 )/iV) - ^(Q 2 )^^^ + O^M)). (6.21) 

Expanding as in (4.4) one finds the saddle point solution 

a NLOSP (Q 2 ,N) 



7 s ™ (a s (Q 2 )/Ar,AT)(-( x ' iVLO ( 7 ^LO ( ^ ( Q2 )/Ar))7V))) | 

exp^ 7 SPArLO (a s (g 2 )/A0dW 

(6.22) 

Further corrections can be calculated as in §4. However, this expression contains some interesting 
information - the dominant contribution to the running coupling corrections to the conformal part 
of the NLO effective splitting function. Calculating T(N, Q 2 ) as a power series in a s (Q 2 ) and 
transforming to x space one recovers all the contributions to the splitting functions in §4. One also 
obtains the term — a s (Q 2 )p 1 ' con - f (a s (Q 2 )$,) which is the transformation of —a s {Q 2 )^ 1 (a s {Q 2 )/N) 
(with the /?o-dependent terms extracted), and contributions to the (ioa 2 {Q 2 )p(a s (Q 2 )^) splitting 
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function. This latter term provides the scale appropriate to use in the NLO conformal splitting 
function using the BLM prescription at NLO [36]. This usually gives different choices for the LO 
and NLO scales, which could be particularly important in this case where the scale choice is so 
important. 

Calculating T(N, Q 2 ) from (6.22) the NLO conformal contribution 

-*s{Q )7 MQ )/N) = -a s {Q ) _ x , (7 o (MQ2)/A0) ( 6 - 23 ) 

comes from the argument of the exponential term. The leading contribution to the /^-dependent 
correction to this comes from the expansion of 



/ / o. v , Xi(7°(« s (Q 2 )/iV)) W 5 

to order a s (Q 2 ) which in T(N, Q 2 ) leads to the term 



(6.24) 



R 2 ^ 2 ^7°(as(Q 2 )/iV) A^(7°(as(Q 2 )/iV))) 2 Xi(7 (« s (Q 2 )/iV)) A ( , 

) dHa s (Q 2 )) { (-^(7°(« s (Q 2 )/iV)))3 | (^5) 

It is easy to check that all other terms of 0(f3oa 2 (Q 2 )) are less divergent as N — > X(Q 2 ) than this 
one, including the contributions due to the /?o-dependent term appearing explicitly in (6.20), which 
are very subleading. Similarly, the contributions from the unknown X2(t) will be very subleading 
unless X2(t) is rather singular at 7 = \. Taking the ln(Q 2 )-derivative of (6.23) and transforming 
this and (6.25) to x space one may find the scale for the NLO splitting function in the same way 
that the scale for the LO splitting function was found in §4. However, comparing (6.23) and (6.25) 
with the terms in (4.7) one notices some similarities. These are not accidental, and a careful 
analysis following the lines of (4.13) to (4.20) leads to exactly the same result as at LO - the scale 
appropriate to the NLO conformal splitting function is given by 



1 

ln(QVA 2 ) - ln(QVA 2 ) = ln(Q 2 /A 2 ) + - (« s (Q 2 )g) 2 . (6.26) 



This exact equality was not at all guaranteed and is a remarkable result, implying the universality 
of this scale choice at all orders. It is also renormalization scheme independent, like the asymptotic 
form of the LO scale choice. It is undoubtedly true that the LO scale and the NLO scale will 
differ for finite x, this depending on the unknown NNLO kernel, but it shows that the asymptotic 
results are very simple and perturbation theory ought to be particularly convergent at small x. 
The NLO scale also matches well with the qualitative predictions obtained from consideration of 
the singularity structure of the full solution, as we will see below. 
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Using this scale at NLO we can investigate the precise effects of the NLO corrections not 
associated with the running coupling, the so-called conformal contributions. To begin with I simply 
remove the /^-dependent terms from (6.13) obtaining 

xP(x, Q 2 ) = a s exp(A£) ( - a s ( -^A ) ■ (6.27) 

Therefore, considering constant for the moment, we see that the NLO correction is both 

numerically large, and enhanced by a power of a 2 ^ compared to the LO. This latter point is 
really expected. Consider a leading order result of the form exp(Aa s ^). When we go to NLO 
the coupling constant a s becomes a renormalization scheme dependent quantity, uncertain by 
0(a 2 s ). In order to be consistent with the renormalization group and produce a result which is 
independent of renormalization scheme up to higher orders the form of the full solution must be 
exp((^4a s + Ba 2 + • • ■)£), where B is scheme dependent. Expanding this about the LO solution 
we get exp(^4a s £)(l + Ba 2 ^ + •••), i.e. the NLO correction is indeed a power of a 2 ^ times the 
LO result, exactly what we see in (6.27). From this argument it is clear that the NLO correction 
should be exponentiated, and we obtain 

xP(x, Q 2 ) = a s .°-°™ exp(A£(l - 6.5a s )), (6.28) 

i.e. we obtain (slightly altered due to the removal of the /^-dependent term) the publicized correc- 
tion to the powerlike behaviour. 

However, we know that a s is not a constant, but runs according to our scale choice at both LO 
and NLO. Indeed, the renormalization group argument above shows that the NLO terms in (6.12) 
which behave like — 0.69a s (3oa s exp(A£) are not of the form we would naturally expect for the NLO 
corrections, i.e. are not just a power of a s higher, do not represent the order of renormalization 
scheme uncertainty, and are not really subleading. Resumming by absorbing them into the definition 
of a s seems the only sensible thing to do. Doing this and using the scale choice (4.21) in the small 
x limit in the expression (6.28) in both the LO and NLO parts gives 

xP(x, Q 2 ) oc j-i exp(1.14(£/a s (Q 2 ))5 _ 3.0/a s (Q 2 )). (6.29) 

MQ 2 )0 2 (£M(Q 2 )) 3/4 

Therefore, it is only the LO part which gives the x dependence in this limit. The NLO part gives a 
Q 2 -dependent normalization change, which can admittedly be large (though using the x — > limit 
of (4.21) tends to exaggerate the size of this at finite x), as expected from the singularity structure 
of the solution of the full NLO BFKL equation. Hence, using this scale choice the log of xP(x, Q 2 ) 

25 



is very insensitive to NLO corrections at small x, and we would expect the NNLO corrections to 
-> as x -» 0. 18 

Therefore, I conclude that the remaining NLO corrections, after running coupling effects have 
been absorbed into the LO expression, are made far less significant by the effective scale used, which 
has been shown to be the same for LO and NLO. However, they are still potentially important at 
small x. As far as comparison with experiment is concerned the interesting question is whether these 
NLO corrections are significant within the current range of data available. In order to answer this 
question it is probably better to adopt a more sophisticated procedure, and look at the evolution 
not of some hypothetical gluon structure function, but of the true physical structure functions. 

7. Small x Structure Functions. 

The previous sections have all considered the calculation of the gluon structure function ob- 
tained by integrating the solution of the BFKL equation up to the virtuality Q 2 . Of course, this 
gluon structure function is not a real physical quantity, though it does, as we shall see, contain 
most of the essential information for physical quantities for asymptotically small x. However, we 
would like to see precisely how the results in the previous sections apply to real physical quantities, 
and how universal they are. 

The generalization of the previous results to real physical scattering processes is quite straight- 
forward. Instead of integrating the upper leg of the gluon Green's function from zero up to Q 2 
we perform the convolution of the Green's function with the scattering cross-section for a probe of 
virtuality Q 2 with a gluon of virtuality k 2 [25]. i.e. (2.6) is replaced by 

r°° 

F l (Q 2 ,N) = J —a hg (k 2 /Q 2 ,a s ( f i 2 ))f(N,k 2 ,Q 2 )g B (N,Q 2 ). (7.1) 

Currently the relevent (Ji tg (k 2 /Q 2 ,a s (n 2 )) are known at lowest nontrivial order for a number of 
quantities. This is order a s for Fl(x, Q 2 ) for both massive [37] and massless quarks [38], F2(x, Q 2 ) 
for massive quarks [25] and (dF2(x, Q 2 )/d\nQ 2 ) for massless quarks [38]. For massless quarks 
the lowest order result for F2(x,Q 2 ) is zeroth order in a s and is infrared divergent, representing 
the unknown nonperturbative quark distribution function. None of the cross-sections are known 

This result for the splitting functions as x — > is x P(x,Q 2 ) ~exp{A LO (Z/a s {Q 2 ))2 - B NLO /a s (Q 2 )) 
where Alo is renormalization scheme independent, Bnlo is scheme dependent and higher order corrections 
are claimed to be negligible. The apparent scheme dependence can be eliminated by noting that the 
leading order result assumed ln(Q 2 /A 2 ) = 3.63(a s (Q 2 )0^ as a; — > 0. Including the full ln(Q 2 /A 2 ) = 
ln(Q 2 /A 2 ) + £>LO+3.63(a s ((3 2 )£)2 , where Blo is renormalization scheme dependent, leads to xP(x,Q 2 ) ~ 
exp(A L o(^/a s (Q 2 ))2 +C L o/{a 2 (Q 2 )) + (B LO - B NLO )/a s (Q 2 )) where C L o and B LO -B NLO are scheme 
-independent. 
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beyond leading order, but all diagrams accounting for the running coupling corrections at NLO for 
the structure functions are contained within the NLO BFKL equation. 

Taking the Mellin transformation of (7.1) with respect to (Q 2 /A 2 ) leads to the simple expression 

F i ( 1 ,N) = a s h hg ( 1 )g( 1 ,N), (7.2) 

where as before £(7, N) = 7(7, N, Q^)gs{N, Q0V7, and /ii, 9 (7) is a function of 7 which is finite at 
7 = and 7 = 1/2. Using the appropriate expression for ^(7, N) the inverse Mellin transformation 
may be performed in the same manner as before in order to give the moment space structure 
functions - considering the running coupling constant BFKL equation, either LO or NLO, expanding 
about the same saddle points leads to an ordered solution in a s . Let us examine the simple case 
of .^(iV, Q 2 ) with massless quarks only. As with the gluon structure function it is impossible to 
actually predict this function due to the unknown (N, Qq) and due to the need to regularize the 
BFKL equation when using the running coupling. However, the previous leading twist factorization 
into an incalculable Qo-dependent function and a calculable Q 2 -dependent function also applies in 
the same way. The function /i_l, 9 (7) is entirely associated with the latter and does not alter the 
previous properties for the case of the gluon - the Q 2 -dependent function is a finite unambiguous 
quantity with a Mellin transformation having leading singularities at 7 = and N = 0. 

Hence, as in the case of the gluon structure function the entirely perturbative calculable quan- 
tity to consider is 

r (o 2 m d Hf L (N,Q 2 )) 

T ^l(Q,N)- ^-^ . (7.3) 

This can be calculated for the case of the running coupling and the LO BFKL equation as in §4, 
with all general results being the same as in this previous case, i.e. one obtains an oscillating series 
in Poa s (Q 2 ) and a very similar apparent scale choice, as we will see below. The changes brought 
about by using the NLO BFKL equation with running coupling are also much the same as when 
considering the gluon. As stated, to get a full solution one should use the NLO BFKL equation in 
the form (6.15). Being instead entirely systematic one may use (6.4), and examine the results only 
up to NLO in q s ( / u 2 ). Doing this one calculates the analogues of (6.10) and (6.9). The latter is 
unchanged while the former is altered by the presence of ^l, s (7) into a different coefficient function 
C^ LO (a s /N, Q 2 /y?). The evaluation of this complete NLO coefficient function is not yet possible 
due to the absence of the NLO correction to (JL,g{k 2 /Q 2 , a s (fi 2 )). However, in order to calculate 
the NLO physical anomalous dimension T^°(N,Q 2 / fi 2 ), the analogue of (6.11), one needs only 
the NLO part of C^ LO (a s /N,Q 2 //i 2 ) containing ln(Q 2 /^ 2 ), which is really provided by the LO 
expression via the renormalization group. 19 Explicitly one obtains 

(7.4) 



19 Equivalently one can use the formulae for the physical anomalous dimensions describing the evolution 
of structure functions in terms of themselves, rather than unphysical partons and coefficient functions, 
given in [37]. 
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Hence, the conformal part of Tll(N,Q 2 / p 2 ) is identical to that of T(N, Q 2 /p 2 ), but there is a 
modification to the term determining the scale. In fact, the additional term, h' L g {j ), is a constant 
at 7° = |, and as such it only contributes insignificantly as x — > 0: the asymptotic scale is 

dominated by d f^ a ( ~^J^^ \ an< ^ ^ s identical to the choice already presented for the gluon 
structure function. h' L s (7°) is important at moderate x, however. 

Taking the transformation of (7.4) back to x-space and eliminating the /^-dependent part 
(including the terms in 7 1 ) by setting the scale leads to a precise definition of the effective coupling 
constant to be used for the evolution of Fl(x, Q 2 ) within this expansion scheme. This is presented 
as a function of x for two choices of Q 2 in fig. 6, and can be compared with the effective coupling for 
the gluon structure function (without the /^-dependent terms in 7 1 ) in fig. 3. Clearly in both cases 
the effect of the change in scale is to reduce the small x coupling, and the effect becomes more 
important as Q 2 decreases and the size of p 1 (a s (Q 2 )^) becomes larger relative to p°(a s (Q 2 )^). 
However, for Fl(x,Q 2 ) the effective coupling at x = is larger than a s (Q 2 ). This is mainly 
due to the — (5/3)/?ox(7) term in 7 1 , but is also influenced by the first nontrivial term in the 
series expansion of /il, 9 (7°) in powers of (a s /N) which is negative. As x decreases the effective 
coupling quickly decreases also, and soon falls below that in fig. 3. This latter point is due to the 
^Po(x 2 (l) + x'{l)) term in 7 1 and the remainder of fiL,g("f) which both act to increase p l LL (a s £), 
and hence increase the scale for the coupling. At x = 10~ 5 the effective coupling for Fl(x,Q 2 ) is 
noticeably lower than that for the gluon, but as x decreases even further the effect of the additional 
terms becomes less and less important, and the couplings converge. 

One can now be rather quantitative about the phenomenological effects of the NLO BFKL 
equation and the choice of scale. Let us first make the simple scale choice (i 2 = Q 2 . In this case 
we may write the physical splitting function as 

(x/a s (Q 2 ))P LL (a s (Q 2 )Z) = p° L L (a s (Q 2 )£) - /WQ 2 >I'f (a s (Q 2 )0 - a s {Q 2 )pY L onS {a s {Q 2 )0 

= Pll(*s(Q 2 )0 ~ a s (Q 2 )p L i ot {a s {Q 2 )^), 

(7.5) 

where each of the p LL (a s (Q 2 )^) may be written as a power series of the form 20 

pIl(*s(Q 2 )0 = f j a n {6ls{Q y r . (7.6) 
n - 

The coefficients for the power series of the various terms in (7.5) are shown in table 1. As one can 
see the coefficients for all the p 1 LL (a s (Q 2 )^) are generally much larger than those for p\ L (a s (Q 2 )^). 

Using the conventional choice of scale then at leading order one would obtain the value of 
(dF L (x, Q 2 )/dlnQ 2 ) by convoluting the first term on the right of (7.5) with F L (x,Q 2 ) itself. As 



20 Actually p^ ™^ has an additional term oc 8(1 — x) / (a s (Q 2 )) which appears in the normal one-loop 
physical structure function. 
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an appropriate choice of F L (x, Q 2 ) at a value of (Q 2 /A 2 ) = 8 (Q 2 ~ 40GeV 2 ) I choose F L (x, Q 2 ) = 
(x/0.1) _0 ' 3 O(0.1 — x). This is a function with the approximate shape of F2(x,Q 2 ) at this Q 2 and 
the 6-function is chosen as a crude model for the approximate (1 — x) 6 fall-off at large x. The result 
for the evolution of Fl(x, Q 2 ) is shown in the upper of fig. 7. It increases very quickly at small x due 
to both the shape of Fl(x, Q 2 ) and the large splitting function at small x. Using the conventional 
scale choice one would then find the NLO evolution by using the whole of (7.5). The effect of adding 
in this very large negative contribution to the physical splitting function is also shown in the upper 
of fig. 7. As one can see the effects are dramatic, largely killing the evolution for x > 0.0001 and 
turning it sharply negative below this. 21 Indeed, the NLO correction is nearly as large as the LO 
result for x ~ 0.001, and becomes dominant as x decreases below this: the perturbative solution is 
not at all stable. Also, although we do not have measurements of Fl(x, Q 2 ) in this range of x and 
Q 2 , similar behaviour would feed through to F 2 (x, Q 2 ), and the NLO prediction is dramatically at 
odds with the experimental data. This is therefore a real physical example of the problems induced 
by the NLO BFKL equation, and is completely independent of factorization schemes and hence 
totally unambiguous (which is not the case for discussions of behaviour of the gluon distribution in a 
given factorization scheme). As we go to lower Q 2 the coupling becomes stronger and the expected 
shape of the structure function becomes flatter. Both lead to the NLO corrections becoming even 
more important relative to the LO, and at Q 2 ~ 10GeV 2 the NLO correction is larger than the LO 
for essentially all x. So we see that the conventional choice for the scale leads to disastrous results. 

Let us consider instead the BLM scale choice for Pll(x,Q 2 ). Absorbing p 1 £^(a a (Q 2 )^) into 
the definition of the scale changes (7.5) to 

{x/a s {Q 2 ))P LL (a s (Q 2 )i) = p° LL (a s (Q 2 )0 ~ (*s(Q 2 )S), (7.7) 

where, as I have already noted, the LO scale is only guaranteed to be exactly the same as that to 
use at NLO as x — ► 0. (Using (6.15) it is easy to show that this is true for Fl(x, Q 2 ) in the same 
way as for the gluon - hL, g {"i) only introduces subleading effects as in (7.4).) The result of the 
evolution using the LO splitting function is shown in the lower of fig. 7. It is a little smaller at the 
lowest values of x than for Q 2 = [i 2 , but only by ~ 15%. This is because until we get to extremely 
small x the LO evolution is largely driven by the first term in the power series of Pll(<5 s (Q 2 )£) due 
to the vanishing of the second, third and fifth terms, and relatively small fourth and sixth terms. 
Hence, the decrease of the coupling is only felt as a single power (and indeed there is an increase of 
the coupling for the highest values of x). The discrepancy between the LO results will increase at 
lower values of x. It will also increase as Q 2 gets smaller and/or as the structure function becomes 

21 Similar behaviour was found for the gluon in a particular factorization scheme (MS) using an incom- 
plete calculation of the NLO anomalous dimension [39]. Using the complete anomalous dimension does not 
alter the qualitative results. 
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less steep. It is when we include the NLO corrections that the more dramatic result is seen. The 

size of these now decreases for two reasons: much of the NLO correction has vanished, having been 

absorbed into the definition of the scale 22 , and the effective coupling is now much smaller at small 

x. The result of including the NLO corrections is seen in the lower of fig. 7. It is now a significant, 

but by no means overwhelming effect. As argued in the previous section renormalization scheme 

consistency implies that these NLO effects should really be exponentiated. The result of such an 

exponentiation is also shown in the lower of fig. 7. It is clearly not dramatic, but does help the 

convergence of the perturbative calculation. The exponentiation will become more important as 

x — > 0. Now that I use the BLM scale choice the coupling at small x is far less sensitive to Q 2 

than for fj 2 = Q 2 and the relative importance of the NLO corrections increases far less quickly as 

Q 2 decreases. As shown for the case of the gluon, at asymptotically small x the effective splitting 

l 

function will behave like exp(1.14(£/a s (Q 2 )) 2 ) and the exponentiated NLO corrections will lead 
to an x-independent multiplicative factor. This factor is potentially quite large, however, and the 
NLO effects must ultimately be treated to obtain the correct quantitative results. Nevertheless, it 
appears as though the LO calculation with the correct scale setting may be quite accurate in the 
current range of x and Q 2 probed by experiment. 

These results regarding Fl(x, Q 2 ) seem very pleasing. However, phenomenologically F 2 (x, Q 2 ) 
is far more important since this is the quantity for which we have a great deal of data [40] [41]. 
One can calculate (dF 2 (x, Q 2 )/dlnQ 2 ) in exactly the same way as Fl(x,Q 2 ) simply by using 
the relevent a 2 ^ g (k 2 /Q 2 ,a s (fi 2 )) which leads to h 2tg {^) in Mellin space. Hence, in this case one 
obtains a direct expression for the evolution of the structure function with respect to Q 2 , rather 
than for the structure function itself. However, inverting the Mellin transformation it is easy 
to see that the expression for (d J r 2 (N,Q 2 )/d In Q 2 ) is identical to that for ,Q 2 ) up to the 

/ii, 5 (7)-dependent effective coefficient function (or in this case anomalous dimension). Whereas 
T L (N,Q 2 ) has a factor of C^ LO {a s /N,Q 2 / n 2 ) up to NLO, (dT 2 (N, Q 2 )/d\nQ 2 ) has a factor of 
7^ LO (a s /iV, Q 2 //j, 2 ). Hence, we can write the physical evolution equation 

^^f 1 =T 2L (Q 2 ,N)F L (N,Q 2 ), (7.8) 

where r 2L (iV,Q 2 ) = ^ LO (a s /N, Q 2 //x 2 )/Cf LO {a s /N, Q 2 /fi 2 ) 23 In this expression all the un- 
known nonperturbative physics associated with (dJ 7 2 (N,Q 2 )/dlnQ 2 ) and J r L (N,Q 2 ) cancels out 
to leave us an entirely perturbatively calculable physical anomalous dimension depending only on 
Q 2 , N and at finite order our choice of \i. 



22 For lowish order in the power series the coefficients for p^(a s (Q 2 )£) and p^ on f (a s (Q 2 )£) are similar, 
but the former begin to dominate at higher orders, i.e. lower x, and become totally dominant as n — > oo 
(x — > 0) as demonstrated by the asymptotic results in the last section. 

23 Again one can use the rules for finding physical anomalous dimensions in [37]. 
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As with C^ LO (a s /N, Q 2 / ^j?) we do not know the NLO off-shell cross-section and hence cannot 
fully calculate 7 2 ¥LO (a s /iV, Q 2 /^ 2 ). Hence, we cannot calculate r^ LO (iV, Q 2 ) fully. However, we 
do know all the the effects at NLO due to the running of the coupling for both C^ LO (a s /N, Q 2 / fi 2 ) 
and ^ LO (a s /N,Q 2 /fi 2 ) and can calculate the NLO contribution to T^°(N,Q 2 ) due to the 
running of the coupling and hence find the appropriate scale to use in the LO expression. This is a 
straightforward, though rather lengthy calculation using the NLO BFKL equation in the form (6.4), 
and expanding the Mellin-space solutions about the saddle point for both (dT 2 {N, Q 2 )/dhiQ 2 ) and 
T L (N,Q 2 ) in order to find the relevent parts of C^ LO (a s /N, Q 2 //j> 2 ) and ^ LO (a s /N,Q 2 / > 2 ). It 
results in the relatively simple expression 24 



^ ^7° /Y -X"(7°) LV giy^iM) ) + 1^(7°) ^ !9 (7>2, 9 (7°) 



'ain(a s )VV2x / (7 ) 7°7 V dl J ,o ' 2 \h Lt9 ( 7o ) h% g { l0 ) 

(7.9) 

where h 2 , g (j) and /^(t) can be found in [38], and (/i 2 , s ( 7 ) /h L , g (l)) = f ^(7) = (3/2) 7 +(l- 7 )- 1 . 
As usual we can take the transformation back to x-space. Using the naive scale Q 2 = /u 2 we obtain 

(x/a s (Q 2 ))P 2L (a s (Q 2 ),x) = 5(1 - x)/a s {Q 2 ) + p Q 2L (a s {Q 2 )0 - ^a s (Q 2 )p\f (« S (Q 2 )£), (7.10) 

where the P 2 l(&s(Q 2 )0 are power series of the form (7.6), and the LO physical splitting function 
has a zeroth order term proportional to a 5 function. As is now standard we can find the correct 
scale by eliminating all /^-dependent NLO terms. This is a little more involved than the previous 
cases, but in the asymptotic limit reduces to exactly the same result. As 7 — ► 2 if we keep only 
the most divergent part in the third term on the right in (7.9) then we have the condition that the 
x-space version of 

dln(a s )\ d-f J j0 dln(a s ) \ 2x'(t) J V d l / 7 o 

(7.11) 

must vanish. Since in this limit ^(^2, 9(7^/^,9(7)) ^ ^5/2 this is precisely the same condition 

as we found for the gluon and for Tl (N, Q 2 ) , and we obtain exactly the same asymptotic scale 
(4.21). Indeed, if we attach any physical process to the top of the gluon ladder we will always 
obtain solutions for physical quantities in the same manner: the physical anomalous dimension or 



24 For reasons of simplicity I have previously defined T2l{N, Q 2 /^ 2 ) with an additional factor of a s [42]. 
This leads to no differences when calculating physical quantities. 
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coefficient function will be determined from the part of the solution which has factorized, is (In- 
dependent, and is influenced by the diffusion into the ultraviolet. Hence, we would always expect 
physical quantities to be controlled by the same asymptotic scale. 

Being more precise we may find the x-space version of (7.9) as a power series in a s (Q 2 )^. 
The coefficient functions for the known p 2L (a s S t ) are shown in table 2. Using these series we can 
solve exactly for the scale down to some finite value of x. The effective coupling to be used when 
calculating the small x evolution of F 2 (x, Q 2 ) in terms of Fl(x, Q 2 ) is actually very similar to that 
for the evolution of Fl(x, Q 2 ) over the whole range of x. They become identical as x — > 0, but are 
only slightly different even as x — > 1. 

If we examine the value of (dF 2 (x, Q 2 )/dlnQ 2 ) for the given input for Fl(x, Q 2 ) using the LO 
physical anomalous dimension we find that the decrease in going from the choice Q 2 = fi 2 to the 
effective scale is a little larger than when examining (8Fl(x, Q 2 )/d\a.Q 2 ). This is simply because 
the terms in the power series for p 2L (a s Q are not as small as those for p° LL {a s £), and so higher 
terms in the series, where powers of the coupling are used, are proportionally more important. Since 
we do not actually know the value of p\^ on ^ {a. s ^) it is impossible to evaluate the NLO effects, with 
or without the scale setting, but I imagine they are of similar importance to the those for Fl(x,Q 2 ). 
They will certainly lead to the same general result, i.e. the LO expression ~ exp(1.14(£)/a s ((5 2 ))2 
as x — > with the exponentiated NLO corrections leading to an x-independent multiplicative factor. 

I note that within this picture there is no way of predicting inputs for structure functions (or 
partons) at some fixed Q 2 . However, since the evolution generates no true powerlike behaviour 
there may well be no growth at x — ► stronger than the soft pomeron. I see no reason to believe 
the values for the intercepts calculated by putting some infrared cut-off on the BFKL equation 
for running coupling, which are both cut-off method and scale dependent. However, at the sort 
of values of x we consider in practice, x = 10 -2 — > 10 -5 , the perturbative evolution can generate 
a rise at small x which appears to be like an effective power over this restricted range in x. In 
broad terms this will not be dissimilar to that generated by the NLO in a s (Q 2 ) evolution, but will 
be different in detail. Perhaps the best method for attempting to predict the shape of a structure 
function at a given input scale is to demand that the general form of the structure functions are as 
insensitive to changes in starting scale as possible [42]. In this way the inputs are determined largely 
by the form of the evolution, and hence the effective physical splitting functions. Since the small 
x evolutions of F 2 (x,Q 2 ) and Fl(x,Q 2 ) are related in a calculable manner this imposes a precise 
consistency requirement on the small x inputs of the two. A more detailed study of study of this 
would be interesting, though an obvious conclusion is that the shape of F 2 (x,Q 2 ) and Fl(x,Q 2 ) 
with x should be roughly the same at all Q 2 and hence at Q 2 (see below). 
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8. Phenomenological Consequences. 

Armed with the small x scale choices for the physical structure functions, it is now possible to do 
a phenomenological analysis. The inclusion of the input singlet quark distribution, or equivalently 
the inclusion of T^iN, Q 2 ) and 1^2 (-^V, Q 2 ) is easy since at LO these are related in a simple manner 
to Tll(N, Q 2 ) and T2l(N, Q 2 ) respectively [37]. Furthermore, they have only a small effect. Much 
more important is the treatment of the LO terms in the physical splitting functions which are less 
singular than 1/x as x — > 0. As shown in [42] a full LO analysis should include all such terms 
at lowest order in a s as well as all terms in the LO small x expansions considered so far in this 
paper. A correct extension of [42], which used the simple scale choice Q 2 = /i 2 , would involve the 
full LO, in a s as well as ln(l/x), physical splitting functions with the scale choice determined not 
only by the NLO running coupling effects considered in this paper, but also by the x-finite NLO in 
a s running coupling effects. 

Consideration of the NLO in a s running coupling effects leads to additional important scale 
changes away from Q 2 = at high x. The evolution of the nonsinglet structure function 
Fj vs (x, Q 2 ) was considered in [43] where it was found that the appropriate scale to use is 

Q 2 = Q 2 ^U(x), (8.1) 
x z 

where k(x) is a relatively smooth function of x from — > 1, k(x) s=s 0.15. Careful consideration 
shows that such a scale change (with some regularization as x — > 1) must be implemented at high 
x for quark driven processes, leading to a larger coupling and quicker evolution. There are also 
nontrivial high x effects in the gluon driven processes due to the NLO in a s running coupling terms. 
This changes the detailed form of the effective coupling already presented in fig. 6 for values of x 
above approximately x = 0.05. For values of x below this the finite x effects on the scale fall away 
quickly. 

One particular consequence of including the full 0{a s ) effective splitting functions is that like 
P gg (z, Q 2 ), Pll{z, Q 2 ) actually leads to a fall with Q 2 for high values of z, the rise only setting in 
when the small x terms become dominant. Hence, the fact that the effective coupling for Pll(z, Q 2 ) 
is actually large at high z increases this negative contribution, whereas the smaller coupling at small 
z decreases the positive contribution, as we already know. This means that, looking at the complete 
convolution leading to the evolution of Fl(x, Q 2 ), the increased negative contribution at high z leads 
to the full scale-fixed LO evolution being reduced compared to the full Q 2 = LO evolution more 
than the consideration of small x effects only in fig. 7 suggests. Inclusion of the high x terms at 
NLO has precisely the opposite effect: this time the positive contribution to the evolution from high 
z due to the 0{a 2 ) terms is enhanced, as well as the known effect of the negative contribution from 
small z being much reduced in size. Hence, the negative NLO correction at small x is significantly 
reduced compared to that seen in fig. 7. Details will be shown in a future paper [44], but the 
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apparent convergence of the perturbative expansion is considerably better even than that implied 
in the previous section. 

I leave a full discussion of the implementation of a full LO in a s and ln(l/x) (denoted by 
LORSC) global fit using scale setting in physical anomalous dimensions to a future paper. 25 Details 
of such a (slightly approximate) fit have already been briefly reported in [27], and here I report the 
most important consequences. 

1. Compared to the most recent NLO in a s (Q 2 ) global fit [47] the quality of the x 2 is improved 
from 1511 to 1339 for 1330 structure function data points. (Constraints from non-structure function 
data, e.g. prompt-photon, Drell-Yan etc. at high x are imposed in the same manner for both.) 
A breakdown of the % 2 for each experiment is shown in table. 3. This extremely statistically 
significant improvement is achieved in all regions of x and Q 2 - the scale choice (8.1) helping at 
high x and the resummation of a s £ terms coupled with the scale choice helping at small x. The 
value of the LO coupling is set at a s (M§) = 0.116, where this LO value is unambiguous, contrary 
to the normal case at LO, because the scale choice has been determined unambiguously. The effects 
of varying the coupling remain to be investigated. A standard NLO in a s fit with BLM inspired 
scale fixing has also recently been performed [48] with less impressive results, particularly at small 
x. 

2. Since the procedure for calculating the evolution is very different from the NLO in a s {Q 2 ) ap- 
proach, predictions resulting from the best fit are significantly altered. For example, the additional 
terms in powers of a s £ in p2 L (a s ^) compared to the NLO in a s {Q 2 ) approach more than compen- 
sates for the decrease in the effective coupling at moderate x and Q 2 , leading to a smaller Fl(x, Q 2 ) 
(very similar to that predicted in [42] if Q 2 > 15GeV 2 ) being required to obtain a similar rate of 
evolution for F 2 (x, Q 2 ). Predictions for other processes, e.g. Drell-Yan production, are potentially 
very different in the two approaches. 

3. There is a failure of the NLO in a s (Q 2 ) approach at small x for Q 2 < 2 - 3GeV 2 . This can be 
seen in two ways. If the gluon (and hence Fl(x,Q 2 )) is required to be positive definite down to 
Q 2 < lGeV 2 then the value of (dF 2 {x, Q 2 )/dlnQ 2 ) becomes too large for Q 2 < 2 - 3GeV 2 [49] (a 
plot can be found in [50]), as can be seen by comparing the data with the prediction from a GRV 
type parameterization [51]. Alternatively, the value of (dF 2 (x, Q 2 )/dlnQ 2 ) can be made correct 
down to ~ lGeV 2 , at the expense of having a valencelike gluon distribution, and hence odd shaped 
F L (x, Q 2 ) (see below), at Q 2 = lGeV 2 , and hence negative gluon and F L (x, Q 2 ) below this [47] [50]. 
Each case demonstrates that the NLO in a s (Q 2 ) approach is breaking down at Q 2 ~ 2 — 3GeV 2 

25 It is also necessary to treat the heavy partons in a consistent manner. The way to do this in the 
context of the full LO physical anomalous dimensions with Q 2 — fi 2 was presented briefly in [45], and will 
be presented in more detail in a future paper [46]. 
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at small x. 26 While this might not seem surprizing since there are many potential reasons for this 
failure (higher twist, higher orders and of course ln(l/x) resummations) , it is a problem not shared 
by the full LORSC fit with the correct scale (even though it is a considerably better fit at small x 
than in [47]). Because the small x effective coupling becomes proportionally smaller compared to 
a s (Q 2 ) as we tend to lower Q 2 , and because, as seen in table 1, the coefficients in the expansion 
of p° LL (a s ^) are small, the evolution of Fl(x, Q 2 ) is slowed down at very small x and Q 2 compared 
to the NLO in a s (Q 2 ) approach. Hence, the Fl(x,Q 2 ) predicted by the global fit does not evolve 
backwards into a pathological form at Q 2 = lGeV 2 . This is shown in fig. 8 where I compare the 
predicted F L (x,Q 2 ) with that obtained from the MRST analysis at Q 2 = 1.2GeV 2 . Clearly the 
shape of the LORSC Fl(x,Q 2 ) is not dissimilar to that of F 2 (x,Q 2 ) at the same Q 2 , while the 
MRST Fl(x,Q 2 ) is rather odd, though it looks sensible by about 2GeV 2 . (The rise at very small 
x in the MRST curve is due to the small quark contribution becoming dominant over the large but 
valencelike gluon contribution.) Evolving downwards the MRST Fl(x,Q 2 ) dips down to negative 
values at about lGeV 2 while the LORSC Fl(x,Q 2 ) will clearly be sensible to much lower values 
(this will be investigated in detail in [44])). Since the effective coupling at small x is so small it 
seems reasonable to believe that the full LORSC calculation should really represent the physics 
down to low Q 2 , as it does, whereas even if the NLO in a s (Q 2 ) approach had worked we would not 
have known why. 27 

Hence, all details of the phenomenology of the scale fixed LORSC analysis seem very satis- 
factory, being a distinct improvement on the standard approach and the LORSC analysis with 
Q 2 = fi 2 . As a word of caution, the analysis presented is still a little approximate, and all quanti- 
tative results are likewise approximate. A more careful detailed analysis will appear soon, though 
it would be very surprizing if the same quality fit were not achieved simply by a slight alteration 
of input parameters and hence very slightly different predictions. 



9. Conclusions. 

I have presented a full discussion of the effect of the NLO corrections to the BFKL equation. 
I have shown that if one resums the ln(A; 2 /^ 2 ) terms into a running coupling constant, as must be 

26 I note that despite reports to the contrary an analysis of data using the leading ln(l/x) terms with 
a s (Q 2 ) does not fail in any more dramatic a manner than this. As shown in [45], using the LO physical 
anomalous dimensions to perform the analysis, rather than some factorization scheme which leads to 
extremely ambiguous results at small x, a fit of even better quality than the NLO in a s {Q 2 ) fit can 
be achieved. The only failings are that the pathological behaviour in the predicted Fl(x,Q 2 ) sets in at 
very slightly higher Q 2 , and of course the NLO corrections using this approach appear to be huge. 

27 A recent discussion of the "Caldwell plot" using the LO BFKL equation with running coupling, though 
with very different techniques from those used in this paper, appears in [52] 
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roughly correct, this alters the whole structure of the solution to the BFKL equation. As previously 

pointed out [15] [21], at leading twist it leads to the solution factorizing into a input dependent part 

which requires regularization, i.e. is infrared renormalon contaminated, and a /c 2 -dependent part 

which is well defined. The degree of uncertainty associated with the input part is shown to have 

exactly the behaviour predicted by Mueller [23] [7]. However, this ambiguity affects the input part 

only, not the whole solution. I note that the evolution part as a function of 7 and TV no longer has 

singularities to the right of zero for either 7 or N, a result which has previously been noted [15] [20], 

but seems to have been universally ignored. Hence, this calculable /c 2 -dependent solution has no 

true powerlike behaviour in either k 2 or x - the hard pomeron intercept is zero. These results 

require no assumptions at all. If one takes the running of the coupling in the BFKL equation 

seriously, the input term is indeterminate unless 13/3oas(Qo)£ ^ 1> an d the evolution term is well 

defined and calculable, and has no true powerlike behaviour. This is not difficult to understand in 

a qualitative manner. It has long been known that the typical virtuality of a gluon in the ladder 

1 

representing the BFKL Green's function has a mean of order k 2 , but a deviation of order (a s £)2 

[22]. I have shown that the diffusion of k 2 into the infrared influences only the input dependent 

solution, the strong coupling then leading to infrared renormalons, while the fc 2 -dependent part is 

influenced only by the ultraviolet diffusion. This means that as one goes to smaller and smaller x 

_i 

the appropriate scale becomes larger and larger, the coupling weaker and weaker (like £ 2), and 
the growth from the ln(l/x) terms is sufficiently weakened by the coupling to destroy the powerlike 
behaviour. 

Using the LO BFKL solution with running coupling I have argued that in order to investi- 
gate perturbatively calculable physics one must investigate physical anomalous dimensions [37], or 
splitting functions, which tell one how unambiguous physical quantities evolve in terms of each 
other, and hence are themselves unambiguous, i.e. independent of factorization schemes or scales. 
This is important when using a small x expansion even at low orders due to large factorization 
scheme uncertainties, but is now vital in order to obtain well-defined, perturbatively calculable 
results. While, of course, it is ultimately necessary to use real structure functions F2(x,Q 2 ) and 
Fl(x,Q 2 ), one may for simplicity work with an unphysical, but unambiguously defined gluon struc- 
ture function G(x, Q 2 ). By calculating the solutions for the Q 2 -dependent factors of the structure 
functions about the saddle points, one obtains ordered power series in Poa s (Q 2 ) for the physical 
anomalous dimensions. While these series appear to be very badly convergent, the coefficients 
oscillate in sign, rendering them summable. I hypothesize that one can approximate the whole 
result by using the BLM scale fixing procedure [24] absorbing the NLO /^-dependent term into the 
definition of the scale used in the LO expression. This results in an effective coupling of the form 
l/(/3 (hi(Q 2 /A 2 ) + 3.63(<5 S (Q 2 )£) 2 )) as x — > 0. For different physical variables the moderate x 
couplings are slightly different but the asymptotic form is universal. It is not guaranteed that this 
choice of coupling is really correct. However, the explicit NNLO calculation supports the procedure 



36 



strongly, and it is also consistent with the qualitative features one knows must be associated with 

the full summation, i.e. it smooths out the powerlike growth in x in precisely the correct manner, 

as well as the picture of ultraviolet diffusion. 

Examining the full NLO BFKL equation I find that as far as running coupling is concerned by 

far the dominant effect is produced solely by the ln(/c 2 /^ 2 ) term. All additional NLO /^-dependent 

corrections lead to modifications to the physical splitting functions which are not only numerically 

_ l 

small, but are reduced by a factor of (a s (Q 2 )(,) 2 . This indicates that it is likely that at all orders 

the ln(/c 2 //x 2 ) terms will lead to the dominant small x effects due to running of the coupling. Indeed, 

at NLO the contribution to the physical splitting function from this term is also dominant to the 

1 

conformal corrections by (a s (Q 2 )£) 2 • The latter are of the form expected from a renormalization 
group argument, i.e. a factor of a s (a s ^) up on the LO expression, while the running coupling effect 
is of an unexpected, more leading form, and essentially demands to be resummed. I also proved 
that if one assumes the dominance of the ln(fc 2 //x 2 ) terms the appropriate scale to use at NLO is 
precisely the same as the LO scale as x — > - a result which was by no means guaranteed to be true 
and seems strongly suggestive of the correctness of the approach. It also implies that perturbation 
theory at small x should be particularly convergent. Using this effective scale choice in the coupling 
I find that the remaining, conformal NLO corrections to the physical anomalous dimensions are 
much more under control than for the scale choice Q 2 = [i 2 due to the smallness of the effective 
coupling as x — ► 0. At all x and Q 2 they are subdominant to the LO result, although they can be 
significant, and in the region of x and Q 2 probed at HERA they are numerically quite small. 

An analysis of data using the full LO physical splitting functions containing both leading 
in ln(l/x) terms and all 0(a s ) terms, with scale fixing appropriate to this combined expansion 
scheme, is very successful. It produces a far better fit to data than conventional approaches, and 
also predicts an Fl(x,Q 2 ) of the same shape as F2(x,Q 2 ) down to Q 2 = lGeV 2 , and possibly 
below. In fact, it seems to work perfectly over the whole range of parameter space one might hope. 
The fit to F2(x,Q 2 ) also leads to predictions for other quantities such as Fl(x,Q 2 ) (difficult to 
measure), F^(x, Q 2 ) (not much different from the standard approach) and Drell-Yan production (if 
the necessary BFKL coefficient functions were calculated). 

Since the coupling at small x is weak, seemingly at all orders, one may be optimistic that it is 
possible to use even LO perturbation theory down to very low Q 2 at small x. Indeed, the prediction 
is that the corrections at NNLO and beyond will be insignificant due to the fall of the coupling 
overwhelming all possible enhancement due to small x terms. However, there are still potentially 
important higher twist (A 2 /Q 2 ) contributions. Nevertheless, the weakness of the coupling may 
make one hope that the small x higher twist effects are strongly suppressed, for example a weaker 
coupling would certainly delay the onset of such effects as shadowing [12] rather significantly. Also, 
I note that within the small x expansion there are no infrared renormalons in the calculation of the 
physical anomalous dimensions. Since renormalons lead to ambiguities which must be cancelled by 
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higher twist ambiguities they are normally taken to be estimates of the size of these higher twist 
contributions - indeed, the scale fixing for nonsinglet evolution at high x [43] does imply renormalons 
of the type already calculated [53] . The absence of the renormalons at small x makes the author at 
least optimistic about the smallness of higher twist effects. Some small x higher twist calculations 
have already been performed [54]. However, since the full physical picture at leading twist only 
appears when performing a full resummed ln(l/x) calculation including running coupling effects, 
a true picture of the higher twist contributions may sadly require similar sophistication (if this is 
possible). I certainly feel that any renormalon calculations performed at fixed order in a s may not 
be representative of the true small x higher twist contributions. If the full LO, with resummed 
terms and scale fixing, analysis is indeed successful to very low Q 2 I would regard this as empirical 
evidence, if by no means a proof, of the smallness of higher twist corrections at small x. 

I have commented on other approaches to the NLO BFKL equation throughout this paper. 
There have also recently been alternative attempts to improve the apparent bad convergence of 
the perturbative series which are somewhat orthogonal to the line taken in this paper. In [10] and 
[11] progress is made by finding resummations which improve the convergence of the expansion 
of the kernel, thus implying a sensible, stable pomeron intercept. I have no argument with this 
approach and believe that for single scale processes it is vital for obtaining a stable expansion for 
general values of x. However, I also believe that for structure functions it leads to effects that are 
completely subdominant to those induced by the running of the coupling. If my assumption about 
the running coupling in the kernel being accounted for by the effective x-dependent coupling in 
physical structure functions has any truth in it, it makes resummations of the conformal part of 
the kernel unimportant since the higher orders are so greatly weakened by the reduction in the 
coupling. Hence, while the work in [10] [11] is certainly interesting, I believe it may be unimportant 
for the real physical results, at least as far as structure functions are concerned. 

Also, there has very recently been a proposal to adopt the BLM scale fixing procedure at the 
level of the eigenvalues of the kernel [55]. This is similar, though not identical to the proposal 
for the change in coupling proposed in [56] when the iVj-dependent corrections to the NLO kernel 
were known. It avoids all the running coupling effects I consider in §4, picking up only those 
in 7 1 in §6, i.e. the ^/3o(x 2 (t) + x'il)) an d — (5/3)PoX(l) terms. This leads to a scale change 
ln(Q 2 /A 2 ) — > ln(Q 2 /A 2 ) +A, where A is very small (and negative). However, the NLO contribution 
to the kernel is renormalization scheme dependent, and this result is in MS scheme. By transferring 
to schemes that the authors reasonably argue are more suited to gluon dominated processes, i.e. 
the MOM [57] or T — > ggg [24] schemes, the scale change at 7 = \ becomes Q 2 ~ 120Q 2 , and the 
intercept becomes ~ X(Q 2 )(l— Aa s (Q 2 )). Hence, the large increase in scale and significant reduction 
in the NLO coefficient leads to a sensible NLO intercept of ~ 0.15 which is not too sensitive to Q 2 . I 
believe the eigenvalue of the kernel is an inappropriate place to make the scale choice since, as soon 
as one introduces the running coupling into the BFKL equation, the whole structure changes. The 
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Q 2 -dependent eigenvalue is no longer a real eigenvalue, as it is at strictly LO, and it no longer has a 
direct physical interpretation. This is identical to the statement that the argument of the exponent 
in (6.8) does not in fact truly represent the full evolution of any physical quantity, is by no means a 
true anomalous dimension, and should not be used for setting the scale. In essence the choice in [55] 
misses the most important results generated by solving the BFKL equation with running coupling 
and looking at physical quantities. This is easily seen by the fact that in any renormalization scheme 
the change in scale using the method in [55] is always of the form ln(Q 2 /A 2 ) = ln(<5 2 /A 2 ) + A rs , 
where A rs is a constant depending on the scheme. Using the BLM method for physical quantities, 
as in this paper, always results in \n(Q 2 /A 2 ) = ln(Q 2 /A 2 ) + B rSii + 3.63(a s (Q 2 )£) 2 , where B rS:i 
depends on renormalization scheme and process. Clearly the ^-dependent term is the dominant 
one at small enough x and contains the most important physics contributing to the scale fixing. 
Note that this contribution is also scheme- independent and the same at NLO as at LO, and that 
the choice of renormalization scheme only leads to subleading contributions to the scale at small x. 
Nevertheless, the type of renormalization scheme considered in [55] leads to a value of B rSi i that 
is rather large. This implies that the details of calculations of structure functions in the current 
experimental range may be sensitive to the renormalization scheme chosen. However, when doing 
a full analysis one should use the same scheme for all physical splitting functions, which will be 
influenced by both gluon and quark dominated processes. There are also further changes to the 
scale due to the running coupling effects at 0(a 2 ), which will be scheme dependent, and potentially 
of similar importance to the differences in B rSt i at the relatively high x values where it is relevant. 
A full understanding of the relevance of renormalization scheme changes needs to take these into 
account carefully. 

Hence, to summarize, I believe that the method of solving for physical quantities using the 

BFKL equation with running coupling and full NLO contributions presented in this paper is the 

best way to proceed for the analysis of deep inelastic scattering at small x. Certainly, the conclusion 

that the running coupling serious alters our picture of BFKL physics, destroying predictivity for the 

input and maintaining it, but smoothing out the powerlike behaviour for the calculable evolution, 

seems to be incontrovertible. More controversial is the proposal that the true physics may be well 

_ l 

described by a coupling which falls as x falls like ln(l/x) 2 . This is strongly supported by current 
finite order expansions, the universality between DIS processes and different orders, the diffusion 
picture, and the general features that the full solution must exhibit. However, it may well be 
possible to validate this more strongly, or invalidate it. Also, the discussion in this paper has very 
firmly used the assumption that the lower end of the gluon ladder is fixed at some low scale, as 
is appropriate for deep inelastic scattering. Further investigation is required in order to consider 
different types of process, although I imagine that the qualitative results will be the same. Overall 
normalization will be infrared renormalon contaminated, since even if there are no small scales in 
the problem the diffusion into the infrared will eventually be important for small enough x, while 
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evolution will be calculable but not truly powerlike. If the general results of this paper are correct, 
perturbative calculations at small x will be very reliable and convergent. They would also explain 
why perturbation theory appears, at least qualitatively, to be working at very low scales at small 
x, but also implies that the standard NLO in a s (Q 2 ) approach is not really quantitatively correct 
at small x. More phenomenological work, including calculation of currently unknown coefficient 
functions as power series in a s ln(l/x), would then be important in order to produce truly precise 
calculations for small x physics. 
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Table 1. 



The coefficients in the power series p l LL (a s (Q 2 )£) = a n (a s {Q 2 )^) n /n\ for the various LO 
and NLO contributions to the physical splitting function Pll{x,Q 2 ). 



n 


P°ll 


l.tot 

Pll 


Pll 


l.con f 

PIl 





1.00 


0.23 


-2.00 


1.57 


1 


0.00 


4.38 


4.15 


1.60 


2 


0.00 


15.87 


11.32 


8.29 


3 


2.40 


13.41 


-16.18 


24.25 


4 


0.00 


86.26 


76.03 


35.31 


5 


2.07 


252.92 


167.34 


140.81 


6 


17.34 


323.08 


-81.51 


377.69 


7 


2.01 


1699.65 


1472.42 


713.25 


8 


39.89 


4338.69 


2665.07 


2553.16 


9 


168.75 


7592.65 


1674.16 


6470.97 


10 


69.99 


33409.13 


28319.16 


14435.29 


11 


661.25 


79427.26 


47284.56 


47746.61 


12 


1945.31 


173361.43 


81792.97 


118560.14 


13 


1717.68 


657395.79 


543255.72 


293414.46 


14 


10643.26 


1527235.16 


927749.64 


905642.90 


15 


25266.78 


3833618.50 


23539999.61 


2256438.84 
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Table 2. 



The coefficients in the power series p 2L (a s (Q 2 )£,) = a n("s(Q 2 )£) n / n - f° r the LO and 
/?o-dependent NLO contributions to the physical splitting function P2l(%,Q 2 )- 



n 


P°2L 


i,P 

P2L 





2.50 


-4.00 


1 


1.00 


9.39 


2 


1.00 


36.60 


3 


7.01 


6.27 


4 


5.81 


239.73 


5 


13.40 


687.03 


6 


58.11 


771.35 


7 


64.74 


5281.50 


8 


196.83 


13213.51 


9 


649.89 


24043.80 


10 


930.65 


111578.92 


11 


3034.70 


265509.09 


12 


8527.87 


613964.05 


13 


15046.02 


2311855.03 


14 


48434.53 


5521425.31 


15 


124600.51 


14458201.96 
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Table 3 

Comparison of quality of fits using full leading order (including ln(l/x) terms) renormalization 
scheme consistent expression, with BLM scale setting and the NLO in a s (Q 2 ) fit [47]. The references 
to the data can be found in [47]. 



'FjXti pri tti PTi t, 

A J-/v l_/ KjL. ±±.±L.\—'±.± U 


data 


y 2 






points 


LO(x) 


MRST 


HI F^ p 

LLL 1 2 


221 


149 


164 


ZEUS F 2 ep 


204 


246 


270 


BCDMS F 2 MP 


174 


241 


249 


NMC F£ p 


130 


118 


141 


NMC Fg d 


130 


81 


101 


NMC Fg n /F£ p 


163 


176 


187 


SLAC F£ p 


70 


87 


119 


E665 F£ p 


53 


59 


58 


E665 F% d 


53 


61 


61 


CCFR F% N 


66 


57 


93 


CCFR F£ N 


66 


65 


68 


total 


1330 


1339 


1511 
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Figure Captions 

Fig. 1. Comparison of the full function exp(7 ln(Q 2 /A 2 ) — X(^f) /{(3qN)) to the case where the 
exponent is truncated at C(7 3 ) along the line ^(7) = 5- I choose ln(Q 2 /A 2 ) = 6 and 
1/(P N) = 2. 

Fig. 2. Comparison of the full function exp(7 ln(Q 2 /A 2 ) — X^) / (J3qN)) to the function appearing 
in the saddle point estimate along the line 5R(7) = j°(a s /N). I choose ln(Q 2 /A 2 ) = 6 and 
l/(PoN) = 2, so j°(a s /N) = 0.384. 

Fig. 3. The effective coupling constant for Np = 4 for the gluon structure function as a function 
of x compared to the constant values at the relevant values of ln(Q 2 /A 2 ). 

Fig. 4. Comparison of the exact NNLO splitting function p 2 (a s (Q 2 )^) with the value pre- 
dicted from the choice of scale in the LO expression. Both terms are weighted by 
(exp(A(Q 2 )0)- 1 . 

Fig. 5. Comparison of the exact NLO /^-dependent splitting function p 1 (a s (Q 2 )£) including the 
corrections from 7 1 with the value of p 1 (a s (Q 2 )0 obtained using the assumption in §4. 
Both terms are weighted by (exp(2a s ((3 2 )£)) _1 for ease of comparison. 

Fig. 6. The effective coupling constant for the physical splitting function Pll{%, Q 2 ) for Np = 4 
as a function of x compared to the constant values at the relevant values of ln(Q 2 /A 2 ). 

Fig. 7. The values of (8Fl(x, Q 2 )/ln(Q 2 )) using the resummed physical splitting functions for an 
input of F L (x,Q 2 ) = (x/0.1)-°- 3 9(0.1-x) at ln(Q 2 /A 2 ) = 8 as a function of x. The upper 
figure shows the LO and LO + NLO results for the conventional scale choice Q 2 = fi 2 . 
The lower figure shows the LO, LO + NLO and LO + exponentiated NLO results for the 
x dependent scale choice in this paper. 

Fig. 8. Comparison of the predictions for Fl(x,Q 2 ) at Q 2 = 1.2GeV 2 from the global fit per- 
formed in this paper and the NLO in a s (Q 2 ) fit in [47]. 
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Expansion about y=\ll 
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full solution 

solution truncated at O(^) 
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Expansion about Saddle Point 
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Comparison of NNLO Corrections 
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